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SUMMARY REPORT 

This report summarizes the fallout pattern scaling relationships that were 
developed in the period 1962 to 1964; the report includes the values of the 
scaling equation coefficients that were derived from selected fallout pattern 
data. The meaning of the scalar wind speed multiplier that is used in the 
scaling equations is discussed relative to computer applications of the scal¬ 
ing system and approximate wind speed adjustment factors for use with wind 
speed averages that may be assumed in such applications are provided. The rel¬ 
ative degree of wind shear inherent in the scaling system parameters is also 
discussed in some detail. Basic equations for relating surface density of 
radionuclides and air ionization rates f including consideration of fractionation, 
surface roughness, and instrument response, are given and discussed together 
with the influence of these factors and others on the limiting values of K 
factors that represent the relative amount of the radioactive sources contained 
within the deduced area covered by the fallout patterns. Scaling equations and 
data are also presented for use in estimating, for any location in the fallout 
region, the time of fallout arrival, the time of fallout cessation, the varia¬ 
tion of the exposure rate (i.e. , air ionization rate in roentgens per hour) 
with time during fallout arrival, and the total exposure from the time of fall¬ 
out arrival to selected later times. 
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FOREWORD 

The major content of this report was developed by the author in the 
period 1962 to 1964 and, in draft form, the material has been available to 
computer programmers at the Stanford Research Institute (Menlo Park, Califor¬ 
nia) and at the American Research Corporation (Fullerton, California). The 
computer programs, in turn, have been available to the Office of Civil Defense, 
their contractors, and others. In some cases programming simplifications and 
interpolation schemes have been added to decrease computing time or for other 
reasons. Since such changes can become iterative with respect to departures 
from the original systems, and since without the original scaling functions, 
the program user has no means of checking the program output, it was requested 
that the original scaling equations be reported for record and for computer 
program verification. 

In preparing the report for publication, a few changes in scaling func¬ 
tions were made, mainly on the procedures for estimating the time of arrival 
anH rate of arrival of the fallout from cloud heights. Also, a few out-of- 
date assumptions, statements, and conclusions were deleted or revised. New 
work reported elsewhere (such as that sponsored by the Defense Atomic Support 
Agency, Department of Defense) since 1964 is not discussed nor is reference 
made to such studies, since the results therefrom were not available for the 
analytical results summarized in this report. 

Some question still exists with respect to the appropriate value of a 
wind speed that should be applied in the calculations to conform with the total 
and angular shear effects included in the derived scaling functions. To clarify 
this question, particle displacement calculations have been added, discussed, 
and a first order wind speed adjustment factor is suggested for use when an 
average wind speed is assumed; the details of the discussion on this question 
are given in paragraph 3.5.4. 
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Chapter 3 

DISTRIBUTION OF LOCAL FALLOUT 


3.1 Background 

3.1.1 The Distribution Proceas 

A very simple descriptive statement of the fallout process Is 
that a cloud of particles is formed rapidly as the result of an explosion 
and that this cloud is then dispersed by the wind and by the force of 
gravity acting on the particles to return them to the earth. Most treat¬ 
ments of distribution of fallout assume that the visible volume occupied 
by the nuclear cloud and stem above the point of detonation within a few 

minutes after explosion more or less defines the volume source of the 
1 2 

fallout particles. One treatment, however, considers the particle source 
volume contained within the air volume 3wept through by the rising fireball. 
In i.lther case the source volumes for the particles depend on total yield 
and, if other than surface detonations are considered, on the height or 
depth of burst. The yield-dependent parameters which are used to define 
the particle source geometry include the cloud height, cloud thickness and 
radius, and, occasionally, the stem geometry, and the time dependence of 
these parameters. 

One important additional factor that is usually considered in the 
fallout distribution process is the spatial concentration of the particles 

in the volume; also qualitative considerations have been given to internal 

3 4 

circulations of the particles by several investigators. ’ The discussion 
in Reference 4 on this circulation is summarized in this chapter. 

The fall trajectory of a particle through the atmosphere depends on 

its own properties and on meteorological factors. The various aspects of 

5 2 1 

these factors have been discussed by Schuert, Anderson, and others. 

The major properties that influence a particle's fall rate through the 
atmosphere are its density, diameter or size, and shape. The major meteo¬ 
rological factors are the wind speed and direction and tne air density and 
viscosity as a function of altitude. 
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The two air properties, of course, are dependent on the air pressure 
and temperature, and these, In turn, change with altitude. The wind speed 
and direction arc also highly variable quantities, since each has both 
spatial and time variations. The vertical motions of the air and particle- 
group diffusion can influence the fall trajectory of particles but arc 
usually not taken into account in the study of the fallout distribution 
process. 


It is of*en found that when the observed cloud rates of rise (or the 
stabilized cloud heights) and the particle fall rates are used to compute 
the time of arrival of particles at locations very close to ground :ero 
where fallout from stem altitudes Bhould predominate, the calculated arrival 
times are quite consistently longer than the observed arrival times. Actually, 
the same discrepancy is often observed for cloud fallout at larger distances. 

The consistency of the longer computed arrival time for particles falling 
near ground zero suggests that when the rising cloud takes on a toroidal 
motion, the larger particles are thrown fiom the gas mass and experience down¬ 
ward accelerations for some rather extended period of time. Because the 
calculated fall rates include only accelerations due to gravity, the computed 
time of fall (neglecting downward accelerations) from the height of the cloud 
would always be longer than the true time. Conversely, when the fall rates 
are used in order to estimate the height of origin of a particle from the 
time of its arrival on the ground (including its rise time), the computed 
height of origin i less than the cloud height. 

This interpretation of the above-mentioned observations of particle 
arrival time may be used to describe, in qualitative terms, the process of 
stem fallout. The rising fireball takes on toroidal circulation as it rises 
from the surface of the ground, and this circulation persists through tran¬ 
sition to cloud form until the internal pressures and temperatures of the 
system approach those of the ambient air, thereby establishing a large-scale 
air circulation. Air and soil particles rise from directly below the cloud 
in a narrow visible stem or chimney, and the surrounding air is entrained 
over the whole length of this stem. This rising material flows into the 


2 
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bottom center of the cloud, and the countercurrent air flow, around the 
periphery of the cloud, is downward. The observable effect, upon occasion, 
is that the mass of particles appears to flow out from the top portion of 
the cloud and then downward. As the cloud approaches its maximum height, 
the circulation pattern apparently rapidly disintegrates or breaks up into 
segregated regions of turbulence under Influence of the ambient meteorological 
forces. 

When the toroidal circulation starts, a particle (or liquid drop) in the 
central region of the cloud would, by centrifugal force, be moved to the outer 
periphery of the cloud and then be accelerated downward at speeds greater than 
the particle's normal fall velocity; it would then be at a lower altitude than 
the cloud when its terminal fall velocity is reached. Even if this centrif¬ 
ugal action and movement to the exterior of the rising cloud did not occur for 
the majority of the particles, they could still fall froir. lower altitudes, by 
virtue of the downward circulation around the periphery of the cloud, than 
would be calculated on the basis that gravity-pull alone was overcoming the 
gross rise r<«te of the visible cloud. 

However, even with toroidal motion, the separation of fallout particles 
by size because of gravity forces is still a valid concept. The smaller 
particles will not move outward by centrifugal forces as far as the larger 
ones in the circulation, and they could be swept back upward through the 
cloud as long as the velocity of the rising all is sufficiently large. This 
type of particle source circulation and ejection can be used to explain the 
observed change in radiochemical composition of different size particles 
discussed in Chapter 2, Reference 6. 

The major radiological factors in the fallout distribution process are 
the fission yield and the variation with particle size of the gross radioac¬ 
tivity carried by particles of a given size. The first essentially determines 
the total radioactivity available for distribution on the particles; the second 
involves the distribution of that radioactivity among particles of different 
sizes. 


3 
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3.1.2 Mathematical Representations 

The original attempt to describe and/or predict the end result 
of the fallou^ distribution process — the fallout pattern — was made by 
C. F. Ksanda and coworkers in 1953. 7 The original scaling method was based 
on the work of Laurlno and Poppoff, 8 which described some fallout patterns 
for low-yield devices derived from observed data obtained during Operation 
Jangle in 1951, The original scaling method was intended for predictions 
or estimates of fallout patterns from yields possibly as high as 10 KT. 

9 

In 1955 the method was expanded to include yields in the megaton range, 
without adequate explicit experimental documentation. This method was sub- 
10 

sequently included in ENW ; however, in the latest edition of this document, 
the fallout pattern scaling is revised. In many damage assessment studies 
of fallout effects, a scaling system is to be preferred over a complex 
mathematical model. 

Mathematical models attempt to establish quantitative values for the 
several fallout distribution parameters mentioned above and to compute the 
activity deposited on the ground at various locations, usually through the 
use of electronic computers. The general approach used and the organizations 
and investigators involved in the development and testing of these models up 
to 1957 is described in some detail by Kellogg. Later developments include 

the work by Anderson, 2,11 by Pugh and Galiano, 2 by Callahan et al., and 
14 

by Rapp, to mention a few of the unclassified repoited studies. A general 
comment on the results might be that none of the models agree with each 
other in several details and that none of the models reproduce very accurately 
all of the few data in the yield range of 1 KT to 15 MT that are experimentally 
available. 

The exact causes of the differences among the various models are difficult 
to isolate for at least two reasons: 

1. Each model is differen* from any other in several of its 
assumptions about para eter values or in its manner of 
handling the many variables mathematically. 


4 
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2. Generally, the reports describing the models do not 
Include sufficient detailed information regarding the 
minor assumptions and the methods used in making the 
computations. 

If the input data in the mathematical models were all more reliably estab¬ 
lished experimentally, many of the differences among them would disappear. 
Whether this would produce better agreement with observations would still 
have to be established. 

In mos mathematical models, after selection of the values of the source 
geometry, trajectory, and radiological factors, the computation is carried 
out by dividing the source geometry for each of several particle size ranges 
into horizontal discs of finite thickness. The location on the ground where 
these "particle discs" land, under the influence of stated meteorological 
conditions, is then calculated. All the activity at each of a series of 
coordinate points is then summed according to the number of different discs 
that land at the point and the amount of activity assigned to each disc. 

So short a summary of the work on the mathematical models should not be 
interpreted to mean that the efforts in the development of the mathematical 
models have been small and unfruitful. On the contrary, much has been 
learned about the fallout process through them, and most of the concepts 
employed by many of the mathematical model developments are covered in the 
remainder of this chapter. But to describe all the work and all the details 
of each model currently in use is not considered to be within the scope of 
this discussion. 

3.1.3 Fallout Pattern Features 


Although observed data on fallout patterns f-om land-surface shots 
of various yields are very meager, the processed data give indication in a 
qualitatitive way of a number of persistent characteristics. For example, the 
general shape of the fallout standard intensity contours (in R/hr at 1 hr) 
from shots in which the wind structures were rather simple, resembles a shadow 
of the mushroom cloud and its stem on the ground. 

5 
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Because oi the shortage of reliable data on the fallout patterns from 
land-surface detonations, any systematic method for scaling fallout patterns 
(i.e., methods for interpolating and/or extrapolating data from one weapon 
yield to another) must take full advantage of all such apparently persistent 
qualitative characteristics of the available patterns. In devising methods 
that can convert the qualitative characteristics to quantitative ones, the 
methods must, of course, be capable of at least reproducing the observed data 
used in obtaining the original scaling relationships, which are given as 
functions of weapon yield. 


Some of the apparently persistent characteristics of the fallout patterns 
from surface detonations are: 

1. In the region near ground zero, the intensity gradient in the 
upwind and crosswind directions is very steep. 

2. The high intensities near ground zero appear as an intensity 
ridge (rather than as a circular peak) displaced in the down¬ 
wind direction. 

3. The length of this high intensity ridge appears to be propor¬ 
tional to the width of the lower portion of the stem. 

4. The peak intensity of the ridge increases with yield in the 
1- to 10-KT yield range and decreases in the 100-KT to 10-MT 
yield range. 

5. The best simple empirical relationship for the variation of 
the intensity with upwind and crosswind distance from ground 
zero, from graphical plots of the data, is that the form Ige" 1 **, 
where I 0 is the ridge peak intensity, k is a constant for 

a given yield, and x is the upwind and/or crosswind distance 
from the upwind shoulder of the ridge peak. 

6. The contours downwind from ground zero appear to be parallel 
to the intensity ridge for its entire length. 

7. At distances greater than the length of the ridge, the 
intensity contours directly downwind decrease with distance 
from ground zero. 

8. At some distance downwind (or perhaps even upwind for very 
large yields), the low-valued intensity contours fan out, and 
the intensities directly downwind from ground zero rise sharply 
with distance and then more slowly with distance to a peak 
value. 


6 




URS 702-1 


9. The distance from ground zero to this downwind peak intensity 
increases with weapon yield. 

10. The magnitude of the peak intensity also appears to increase 
continuously with yield. 

11. The distance between the lower valued contours appears to be 
related to the width of the cloud (not considering wind shear 
differences), and the maximum width seems to occur farther 
downwind than the peak intensity. 

12. The variation of the intensity with downwind distance from 
this outer pattern peak can be approximated within reasonable 
limits of error, if the form of the wind shear pattern is not 
complex, by a function of the form I 0 e -mx , where 1 Q is the 
peak intensity, m is a yield-dependent parameter, and x is 
the downwind distance from the peak. 

The above-listed fallout pattern characteristics are based on a combina¬ 
tion of experimental observations and analyses of field test data. The most 
reliable persistencies appear to be those numbered 1, 2, 4, 6, 7, 8, and 9. 

3.1.4 Simplified Fallout Scaling System for Land-Surface Detonations 

The fallout scaling system described here was developed for 
estimating standard intensities, potential exposure doses, and other radio¬ 
logical quantities by use of both manual and machine computational techniques. 
The system is based on corrected experimental data, on empirical relationships 
among the geometrical arrangement of the cloud and stem as the source of fall¬ 
out particles, and on several of the observed features of the fallout pattern 
of radiation intensities on the ground. In the system, the cloud and stem 
dimensions are stylized as simple solid geometric configurations to facilitate 
the use of algebraic relationships among the model parameters and the dependence 
of the parameter values on weapon yield. 

In making estimates of the hazards from fallout, for the purpose of 
establishing the nature and required degree of protection against these hazards, 
two major quantities requiring evaluation are (1) the exposure dose levels 
that can result at different distances from the detonation and (2) the land 
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surface area In which the exposure dose Is greater than a stated amount. To 
make these evaluations requires estimates of the amount of fallout that 
deposits at various locations, the time at which the fallout arrives, and 
the rate of its arrival. 

Such general evaluations of radiological hazard levels and of the pro¬ 
tection requirements for radiological countermeasures in defense planning 
must first consider the possible levels of effect (or hazard) and, in a 
generalized manner, the feasibility of methods for protecting against these 
levels of possible hazard. For these purposes, a rather simplified fallout 
scaling system can serve; no precise or accurate prediction of fallout under 
specified detonation and wind conditions is needed or is possible, even with 
the most complicated fallout models at their present stage of development. 
Therefore in the following discussion the presentation is limited to the 
description of a simplified version of the fallout distribution process. 

The mathematical derivations of the simplified fallout scaling system 
are designed to describe the fall of particles of different size-groups 
from a volume source in the air; the boundaries of that source are assumed 
to depend only on weapon yield. The problem is to describe mathematically 
the dependence of the fallout pattern features, in space and time, on (a) 
the cloud and stem geometry, (b) the particle fall velocities, (c) the wind 
velocity, (d) the radioactivity-particle size distributions, and (e) the 
weapon yield. 

The geometrical configuration of the cloud for the scaling system is 
an oblate spheroid, and the configuration of the stem is a frustum of an 
exponential cone or horn whose larger base is approximately adjacent to the 
bottom of the spheroid. The fall of particles from each of these source 
volumes is considered separately. The mathematical description and detailed 
assumptions used in the development of the model are given in Reference 4; 
only those details n< euc J for use of the derived pattern scaling functions 
are presented in the fallowing discussion. 

8 
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Some of the pattern features of interest along the downwind axis (Y = 0) 
of the Idealized fallout pattern are shown in Figure 3.1 as a schematic inten¬ 
sity profile. The numbers shown in the figure correspond to the numerical 
subscripts of the scaling functions. The evaluated scaling functions ;'cr 
these and other quantities are given in the following paragraphs as summarized 
from the data given in Reference 4. In the scaling system, the fallout pat¬ 
terns for the particles falling from stem altitudes and from cloud altitudes 
are computed separately; at locations where these two patterns overlap, the 
computed standard intensities are then added together. ThiB is illustrated 
in Figure 3.1 by the dotted line between X 3 and X g . The standard or reference 
condition for all the fallout pattern scaling functions presented in the 
following paragraphs and sections is 100 percent fission yield; the radioactive 
components assumed to be present are discussed in Paragraph 3.4.3. 

3.2 Fallout Deposition from Stem Altitudes 
3.2.1 Ground Zero Intensity Ridge 

The ground zero intensity ridge is depicted in Figure 3.1 by the 
downwind distances and and the intensity 1^ The dependence of X 2 
and Xg on the average wind speed and weapon yield is represented by 

X 2 =“ 0.0327W°' 230 jv^ - 3.96W°‘ 128 Jmiles ; W = 30 to 10 5 KT (3.1) 

and 

X 3 =- 0.0327W°’ 230 £v w + 3.96W°’ 128 Jmiles; W = 30 to 10 5 KT (3.2) 

for in miles per hour and W in kilotons total yield. The values of 3 
are estimated from 


* Approximate equalities indicate smoothed logarithmic scaling functions; equal 
signs are used for the basic scaling functions assumed in the derivations. 

9 
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l 2 3 = 3,02 X 1(>6v w W W“° -519 R/hr at 1 hr; W = 30 to 10 5 KT (3.3) 

where I 2 3 is the standard intensity for a 100 percent fission yield weapon, 
and n Is an empirically derived yield-dependent parameter whose value is 
calculated from 


= 0.821 - 0.0191 log W; W = 30 to 9,000 KT 


(3.4) 


- 1,200 - 0.115 log W, W - 9 x 10 3 to 10 S KT (3.5) 

The effective fall velocity vector for the median-diameter particles 
landing at (X 2 + X^l/2 is given by 

v = €8.3W" 0,080 ft/sec; W = 1 to 10 5 KT (3.6) 

2,3 

The apparent height of origin of the particles with the fall velocity vector 
given by Equation 3.6 can be estimated from 

Z 2 3 ^ 804 * lo3w °' l5 ° feet! W = 30 to 10® KT (3.7) 


3.2,2 Downwind Intensity Contours 

The downwind distance to the location of Figure 3.1 on the 
center line of the stem pattern is given by 

X, =" 0.316v W 0-203 miles; W = 30 to 10 5 KT (3.8) 

4 w 


The standard intensity, for a 100 percent fission yield weapon, associated 
with the downwind distance designated by X 4 is given by 


I = 15.0/v 


(3.9) 
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The downwind distances, X , to other Intensity values. I , between I and 
c* * * s’ 2,3 

1^ on the pattern centerline from to X^ are estimated from 


(X - X > log (1 /I ) 


(3.10) 


The effective fall velocity vector of the media-diameter particles 

landing »t X, is given by 
4 

v = 13.8W -0 ' 082 ft/sec; W = 1 to 10 5 KT (3.11) 

f 4 

The apparent height of origin of the particles with the fall velocity 
vector given by Equation 3.11 can be estimated from 

z 4 - 1.57 x 10V- 121 feet; W = 30 to 4 X 10 4 KT (3.12) 


or 


4 o no 4 e 

z 4 ~ 1.62 X 10 W - 11 feet; W = 4 X 10 to 10 KT (3.13) 


3.2.3 Upwind Distance to the 1 R / hr at 1 hr Contour 

The upwind distance to the 1 R/hr at 1 hr contour (i.e., 1^), 
designed by the numeral one in Figure 3.1 is estimated from 

X x = -0.695W 0,319 miles; W = 30 to 10 5 KT (3.14) 


for an average 15-raph wind speed. For other average wind speeds, this distance 
on the pattern center line is calculated from 

0 337 5 

* l = X 2 - 0.174W • log I 2 3 miles; W = 30 to 10 KT (3.15) 

where X g and 1^ 3 are given by Equation 3.1 and Equation 3.3, respectively. 


12 
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The distances to locations between and X g tor the Intensity values, I s> 
other than those ol 1^ and I 2 3 are estimated from 


(X - X ) log 1/1 

X' = X. + -=- rz — -- 

C 1 log 1 jj j/Ij 


(3.16) 


where X^ is the distance from ground zero to the Intensity value I g along the 
center line of the pattern. 

Because the upwind fallout intensity contour locations are based only on 
empirical data, the scaling system cannot be utilized to estimate the diameters 
of the particles landing at locations between X 1 end X 2 . 

3.2.4 Stem Pattern Half-Width 


The (maximum) half-width of the stem fallout pattern, designated 
by Y ., Is the lateral distance from the center line of the stem pattern to 
tnc 1 R/hr at 1 hr contour; It is located at the downwind distance X g given 
by 


(X - X ) > Y 


(3.17) 


or 


X s = X 2 ; (X 2 - X t ) S y g (3.18) 

The value of Y g (15), for an average 15-mph wind speed Is estimated from 

Y g (15) = 0.316W 0,400 miles; W = 1 to 10 5 KT (3.19) 

For other wind speeds, the stem pattern half-width, Y g (v^), is calculated 
from 
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Y (v^) = Y s (15)S s (V w ) (3.20) 

In which S (v ) is a wind velocity correction factor. The values of S (v ) 
s w s w 

are given as a function of weapon yield for several selected average wind 

speeds in Figure 3.2, the value of S (v ) for wind speeds other than those 
given in the figure can be obtained graphically from a linear plot of S^v^) 
against wind speed for a given weapon yield. 

The lateral distances from the pattern center line to other intensity 
values between I 2 g and the 1 R/hr at 1 hr contour (100 percent fission yield) 
at the downwind distance X g are estimated from 


o W loS (I 2,3 /1 s ) 
Y 0 = log l^g/Ij 


(3.21) 


3.3 Fallout Deposition from Cloud Altitudes 
3.3.1 Downwind Intensity Contours 

The downwind intensity contour locations on the center line of the 
fallout pattern produced by particles falling from cloud heights are defined 
in terms of the intensities designated by I 5> I g , Iand I g at the downwind 
distances X,., X g , X ? , and X g , respectively, as shown in Figure 3.1. 

The dimensions and height of the assumed elliptical source volume for the 
particles falling from the cloud starting at times of about 3 to 10 minutes 
after detonation, as derived from the data reported by Pugh and Galiano 1 " and 
by Schuert,^ are represented by the following yield-dependent functions: 

r = 0.464W 0-431 miles; W = 1 to 10 5 KT (3.22) 

b = 0.265W 0 ' 300 miles; W = 1 to 10 5 KT (3.23) 

h = 3.18W°' 164 miles; W = 30 to 10 5 KT (3.24) 
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r 


i 


where r is the cloud radius, b is the cloud half-thickness, and h is the 
altitude of the center of the cloud. To a first approximation, r represents 
the cloud radius at the time the cloud top reaches its maximum height ; the 
visible cloud radius continues to increase for some time thereafter. 

In the fallout pattern scaling system, the variation of the downwind 
distance with weapon yield of the selected contour locations is mainly 
accounted for by setting the various characteristic distances proportional 
to the cloud height for a given average wind speed. For the locations under 
the cloud, some correction to allow for a decreased effective cloud radius 
is included in the scaling function. An additional, but smaller, dependence 
on weapon yield of the distances to the selected contour locations is included 
to account for a gradual shift in the radioactive content and in the average 
diameter of the particles that land at the designated locations. This change 
in radioactive content is apparently caused by a gradual increase in the 
proportion of the radioactive nuclides carried by the smaller particles as 
the weapon yield Increases. 


The scaling functions for estimating the distances and the intensities 
for the selected contour locations are given in two sets, depending on 
whether the downwind distance is less or greater than the cloud radius. For 
average wind speeds less than a given speed, the distances X r and X g are less 
than the cloud radius, r, the scaling functions for the3e distances and wind 
speed limits for a giver, weapon yield are given by 

X 5 = 0.186W°* 254 (v w - 2.49W 0 * 177 ) miles; v^ S 3.28W 0,245 mi/hr; 

5 

W = 30 to 10 KT (3.25) 

and 

X g = 0.179W 0,309 (v w - 2.59W 0,122 ) miles; v^ £ 2.04W 0,231 mi/hr; 

W = 30 to 10 5 KT (3.26) 
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The downwind distances to all selected locations that are greater than 
cloud radius are represented by 

X = v X®W n i miles (3.27) 

i w i 

in which X^^ and n^ are constants, and the subscript i represents one of the 
selected locations; derived values of X A and n^ are summarized in Table 3.1. 

The standard intensity for each of the selected locations is calculated 

from 

1 ± = K°W ra i log cp i (v w ) R/hr at 1 hr (3.28) 

in which and are constants, and cp^fv^) is given by 

Cp 1 (v w > = 

V w 5 v i r/h ! w = 30 to 1q5 KT (3.29) 


<v w + 

VV = - !L — 

<v„ - v t r/h) +^[ V JrV 2 + <* w “ v^/h) 2 

v S v r/h ; W = 30 to 10 5 KT (3.30) 

w i 

in wnich v^ is the average value of the fall vector for the particles deposited 
at the location designated by 1, and v^m) is the minimum value of the fall 
vector of the deposited particles. The value of l g obtained from Equations 
3.27 and 3.28 results in 
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1. = 15.0/v 


(3.31) 


for v In miles per hour. 

w 

The empirically derived values for the constants of the Intensity scaling 
functions for each of the selected contour locations are also summarised in 
Table 3.1. The median diameter of the particles deposited at each of the 
selected locations is estimated from the v^ values calculated from the func¬ 
tions of Table 3.1 (after multiplying by 1.467 to convert the values from 
mi/hr to ft/sec) using data given in Reference 4, along with h the cloud 
center height, as the height of origin of the particle source. 

One of the major characteristics of the fallout patter.i scaling system 
for the fallout from the cloud altitudes is that it specifies a peak in the 
downwind intensity profile at X ? . The intensity, I 7 , gradually increases 
with weapon yield approximately proportional to W 0 * 4 . No experimental data 
are available for testing the reliability of the estimates of I ? obtained 
from extrapolation of the scaling functions to very high weapon yields. 

The downwind distances to contours of other standard intensities on the 
cloud fallout pattern center line are estimated from 


(X - X ) log (I /I c > 

X' = X e + — - --- ■ — —±- ; 1 S I 

c 5 log (Ig/I 5 > s 6 


(3.32) 


(X - X ) log (1 /I > 


log (I/Ig) 


I s I £ I_ 


(3.33) 


cr 


(X - X ) log (I /I 

— Y + ■ — - ■ - - T < T 

7 log (I ? /I 9 ) ’ s 7 


(3.34) 


where X^ is the downwind distance to the Intensity, I g , for the distances less 
than X ? , and X^ is the distance to the intensity, l s> for distances beyond X ? . 
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3.3.2 Maximum Pattern Half-Width 

The maximum pattern half-width for fallout from cloud altitudes 
is designated as Y g , and the downwind distance to Y g is designated as X g 
(see Figure 3.1). The empirically derived scaling function for X g Is given 
by 

X g = 0.325VJ* 0,315 miles; W = 30 to 10 5 KT (3.35) 

The crosswind distances to given contours in the fallout area depend, 
first, on the lateral displacement of the particles during the rise of the 
cloud; second, on the wind directions at all altitudes from the bottom to 
the top of the cloud; and third, on the wind speeds. 

The wind speed has two effects on the lateral displacement of an inten¬ 
sity contour. One is the horizontal displacement of particles with wind speed 
because of the relative horizontal distance traveled In a given period of time. 
The other is the decrease In surface density of a given size group with wind 
speed because of the change in the angle of the particle trajectory. Hence, 
even for the case in which the wind direction is the same at all altitudes, 
a change In wind speed results in a change in the maximum crosswind distance 
of a given intensity contour. 

The values of Y°(15) for the maximum lateral distance from the pattern 
center line to the 1 R/hr at 1 hr contour for an average wind speed of 15 
miles per hour, as derived from both observed data and summations of the 
fraction of the radioactivity contained in the fallout pattern (see Section 
3.4), are summarized in Table 3.2. The computed pattern widths include the 
effect of lateral wind shear contained in the original data; this effect is 
discussed in Paragraph 3.5.4, Approximate scaling functions for Y g (15) are 
as follows: 

Y g (15) =“ 0.518W 0,615 miles; W = 30 to 750 KT (3.36) 
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Table 3.2 

SUMMARY OF CALCULATED Y°(15) VALUES FOR SELECTED WEAPON YIELDS® 


W 

(KT) 

5 

10 

20 

50 

100 

200 

500 


Y s< 15 > 

(miles) (MT) 

1.90 1 

2.28 2 

3.18 5 

5.76 10 

9.10 20 

14.0 50 

23.5 100 


Y°(15) 

(miles) 

33.6 
40.9 
53.4 

64.7 
78.0 

101 

123 


a For 100 percent fission yield 
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and 


Y g (15) =" 4.75W 0-283 miles; W = 750 to 10 5 KT (3.37) 

The variation of Y g with wind speed (for a given wind direction) is deter¬ 
mined relative to Y g (15) for a wind speed of 15 miles per hour. The representa¬ 
tion for the variation is 


W = Y 8 (15)S<V w ) 0.38) 

in which S(v ) is the relative shear factor due to wind speed only. 

The values of S(v ) determined from the fallout scaling system parameters 
for different wind speeds are essentially independent of weapon yield. The 
indicated value of Y g (v w ) is for the particle groups falling at the downwind 
distance, X g ; the associated intensity contour that passes through the loca¬ 
tion of Y g (v^), x g » ls equal to I Q for the same wind speed. The intensity 
at the location is thus 1 R/hr at 1 hr when the wind speed is 15 miles per 
hour. Values of S(v^) at several wind speeds, and the associated intensities, 
are given in Table 3.3; the tabulated values of S(v ) are represented approx¬ 
imately by 


S(v ) *** 0.360(1 + 26.7/v ) ; v = 10 to 22.6 mi/hr (3.39) 

and 

S (v w ) “ 0.426(1 + 19.0/v^); v^ = 22.6 to 50 mi/hr (3.40) 

Combining Equations 3.39 and 3.40 with Equations 3.36 and 3.37 gives, for 

VV' 

Y g (v w ) =“ 0.186W 0-615 (1 + 26.7/v^) miles; v^ = 10 to 22.6 mi/hr (3.41) 

Y g (v w ) =" 0.221W 0-615 (1 + 19.0/v w ) miles; = 22.6 to 50 mi/hr (3.42) 
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Table 3.3 

SUMMARY OF THE RELATIVE WIND SPEED SHEAR FACTOR, S(v ), 
FOR THE FALLOUT PATTERN MAXIMUM HALF-WIDTH AND * 
ASSOCIATED STANDARD IONIZATION RATE 
FOR SEVERAL WIND SPEEDS 


Wind Speed 
(ml/hr) 

10 

IS 

20 

25 

30 

35 

40 


s <V 

1,325 

1.000 

0.840 

0.750 

0.695 

0.658 

0.630 

0.605 


(R/hr at 1 hr) 
1.5 
1.0 
0.75 
0.60 
0.50 
0.428 
0.375 
0.333 
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for W values between 30 and 750 KT; and 

Y g (v w ) “ 1.71W 0,283 (1 + 26.7/v^) miles; v = 10 to 22.6 mi/hr (3.43) 

Yg(v^) =*■ 2.02W 0,283 (1 + 19.0/vO miles; = 22.6 to 50 mi/hr (3.44) 


5 

for W values between 750 and 10 KT. 

The value of Y g for the 1 R/hr at i hr contour for a given wind speed 
(100 percent fission yield) is estimated from 


Y^(l) 


Y g ( v w> log I ? 
log d 7 /I 9 > 


(3.45) 


where I ? and I g are the Intensities at X ? and X 0 (Y = 0) for a given value of 


The maximum lateral distance to other intensity values between I ? and 
the 1 R/hr at 1 hr contour at X , Y°(l) is estimated from 


The downwind distance to Y g (I) is given by 

(X - X ) log (I_/I ) 


log (I/I 9 > 


3.4 Idealized Fallout Pattern Shapes and Areas Within Contours 
3.4.1 Stem Pattern Contours 


With the simplified fallout scaling system and the stylized down¬ 
wind intensity profile, it is convenient to construct contours by using simple 
geometric forms which approximate to some degree the true shapes of the contours 
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and, within reasonable limits, account for the fraction of the activity 
produced that falls back to earth within their areas. 

The assumed shape of the intensity contours for the fallout from stem 
altitudes is a combination of a half-ellipse and a partial circle having its 
center at X 2 , providing the distance (X 2 - X^) is equal to or larger than the 
distance Y for the same contour of the stem pattern. For these conditions, 
the major axis of the half-ellipse on the center line in the downwind direc¬ 
tion is (X^ - x g ) where is equal to X^ + [(X^ - Y^)“ - Y^)'**] 1 ^ 2 and the 
minor axis is Y° the ellipse is centered at x g on the pattern center line. 

The equations for the assumed intensity contour shapes, accordingly, 
are given by 


(X - X 2 ) 2 + Y 2 = <X 2 - X^) 2 ; X = -X^, to x g 


(3.48) 


and 



(3.49) 


The area contained within a given standard intensity contour, for the above 
assumed contour geometry, is given by 


4 - 


sin Y /(X. 


1.5708Y[X c - Xj 


- 0.5708T^(X 2 - X^) 2 - Y 2 ; (X 2 - X^) 2 Y° 


(3.50) 
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or 


log (I 2 3 /I s ) 


[3.1416(X 2 - X/ 

sin _1 Y (v )/(X„ - X.) 

1 _ S W Z 1 

log x 2.3 


| l0g I 2,3 

180 


1.5708Y (v )(X. - X.) 0.5708Y (v 




Y (v )}log (I /I ) 


f 1.5708 Y s (v w )(X 3 



(X 2 - X x ) 2 Y s (v w ) 


(3.51) 


The total activity contained within the stem fallout pattern, obtained 
from the integration of I da (where da is calculated from Equation 3.51 over 
the intensity limits from 1^ 3 t0 I 4 > is B iven by 


*2,3 

[l.l86(X 2 - X x ) 2 | 

sin _1 V s (v w )/(X 2 

log *2,3 

[ iog j 2 ,3 : 

I 1 180 


0.593Y s (v w )(X 4 - X 3 ) 0.2154Y s (v w ) ^(Xg - X^ 2 - 

+ lo * < " 2l 3 /I 4 ) l0g *2.3 


3 


+ 0.682Y s (v w )(X 3 - X 2 )| ; (X g - X^ 2 Y^v^); 1^ » 1 4 (3.02) 

For some conditions of wind speed and weapon yield, the maximum half¬ 
width distance, Y°, will exceed the distance (X 2 - X^). As an example, 

Y (v ) for the 1 R/hr at 1 hr contour exceeds the distance (X - X ) for 
s w 2 1 

weapon yields larger than 22 NTT and a 15 mile per hour wind speed. The 
assumed contour shape for the case where Y^ 2 (Xg - X^_) is that of two half¬ 
ellipses, both centered at the point Xg on the pattern center line. The 
upwind half-ellipse has a major axis given by Y and a minor axis given by 


26 





UKS 702-1 


(X 2 - X^); the downwind half-ellipse has a major axis given oy (X c - X 2 > 
and a minor axis given by Y . The general equations for the idealized 
intensity contour shapes are 


- = 1 ; X £ X„ 


(X- - X')“ (Y ) 


(3.53) 


(x - x 2 ) 2 
«c - V 2 



(3.54) 


The area contained within a given standard Intensity contour, for the 
two half-ellipses, is 

a g = 1.5708Y(X c - X^); Y° * (X g - X^) (3.55) 


or 


1.5708Y s (v w ) log <I 2 .3 /I s ) (P X 2 " V <X 4 " V 

% = io 8l 2 ~ 3 ' j[ lo 8 *2,3 + lo « »2.3 /I 4 ) ^ 

+ (X 3 - X 2 )j; Y s (v w ) 2 (X 2 - X x > (3.56) 

The total activity contained within these contour areas is given by 


0.682Y (v )I 


6 2,3 

' X i> : I, 
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3.4.2 Cloud Pattern Contours 


The assumed shape of the intensity contours for the fallout pattern 
formed by particles falling from cloud altitudes is described by two half¬ 
ellipses smoothly joined at the downwind distance of the maximum width of the 
contour. The locations of the maximum contour widths fall on a line joining 
X ? on the pattern center line (Y = 0) and the point Xg(I),Y g (I) (see Equations 
3.46 and 3.47). 


The general equations for this idealized contour shape are 


[X g (I) - X] 2 y 2 
[x g (D - x ;] 2 y 2 (i) 


1; X <: X g (I> 


and 


[X - X g (I)] 2 
[X c - x 8 (i)] 2 


1; X 2 X g (I> 


(3.59) 


The area contained with a given standard intensity contour is given by 


a = 1.5708Y (X - X') (3.60) 

c See 

in terms of the designated distances; in terms of intensities and scaled 
pattern features, the areas within contours are given by 


1.5708Y 8 (v w )log(I 7 /I s ) f(X 9 " X 7 )log(I 7 /I s ) (X ? - X 6 )log(I s /I 
: = log <i 7 /i 9 > log a 7 /ig) log a 6 /i 5 ) 

* «7 " *5>j ' h 5 >6 
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and 


’.570eY 8 (v v )log(I 7 /I s ) 
'log <I 7 /I 9 > 




- X 7 )log(I 7 /I s ) 
log U 7 /I 9 > 


(X ? - X "->g(I g /I 6 ) 
io 6 CI 7 /I 9 ) 


be 



(3.62) 


The activity Integral for the cloud fallout pattern, designated A , can 

c 

estimated from 


1,364Yg<v^)| (X 8 -X,> (X 7 -X 6 )| I . 

c ■ lex 1 7 |log <I 7 /I 9 ) log <I 7 /I 6 )| j n 7 ' 

0.6822Y e (».)I 6 f (X, - X 7 ) f 

* —nn;— pnvv L 2102 <VV * ° e68 J 

T [ lo * <VV - log <I 6 /I s ) + °- 868 ] + X 7 - X sj 


Ifl) - h <W 


6 6 ‘5 

for the case where I. » I . 


3.4.3 Summary of Pattern Features and Activity Summations 


The fission-product radioactivity produced by a detonation is about 
23 

1.4 X 10 BW fissions where W is the total yield in KT and B is the ratio of 
fission to total yield. Some fraction of this total activity is contained 
within ^ he fallout pattern. If the activity in a pattern is summed over the 
fallout area from the central high intensities down to a stated low-valued 
intensity contour enclosing the largest area, the fraction accounted for 
increases with yield. An ionization rate or intensity sum of a fallout pattern 
made in this way does not account for the activity deposited on the ground at 
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lower Intensities than the selected "lowest” contour, nor does it account for 
the fraction carried away on very small particles as world-wide fallout. The 
sum or integration of the activity over the fallout area can be used, however, 
to determine the fraction of the activity that is accounted for in the fallout 
pattern out to a stated low-level contour; the equations of such a summation 
for the simplified fallout scaling system are given in the two preceding 
sections. 

The amount of activity accounted for in a fallout pattern is defined by 

_ 23 f ° 

C(1)K(1) 1.4 X 10 BW = \ i g da (3.64) 

0 

in which C(l) is the true fraction of the number of fissions accounted for, 
and K(l) is the average value of the ratio of R/hr at 1 hr to the number of 
fissions or kiloton equivalents per unit area for all the activity within 
the area, a . The average or accumulated 1rac.ionation number, for the 
radioactive elements accounted for, ma> be defined through K(l) by 

KO) = 5q Jr fp (l) i fp (l) + i t (l)J (3.65) 

in which 6 is an instrument response factor, q is a terrain attenuation factor, 
r^ p (l) is an average gross fractionation number for the fission products appli¬ 
cable to the ionization rates at H + 1, i. (1) is the H + 1 air ionization 

fP 

rate at 3 ft above an infinite ideal plane for a unit surface density of 

unfractionated fission products spread over the plane, and i^d) is in the 

same units for neutron-induced activities. The units of measure for K(l) 

are the same as for i. (1) and 1,(1) and usually are given in terms of either 
fp i 

R/hr at 1 hr per fission/sq ft or R/hr at 1 hr per KT/sq mi. It may be noted 
that the true air ionization rate is obtained when the value of 6 is set equal 
to one. 
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The empirical constants for the scaling functions of the simplified fall¬ 
out scaling system apply to the reference condition of 100 percent fission 
yield (B = 1), including 0.8 neutron captures per fission by U-238 to give 
the appropriate value of I in R/hr at 1 hr. 

Also, the I g values (or standard intensities) correspond to radiac measure¬ 
ments taken at 3 ft above an extended open area contaminated with fallout as 
it existed when the measurements were taken. The reference radiac instrument 
for the I g values is the AN/PDR-39(T1B) portable radiac; it has a geometric 
and photon energy response at very nearly 75 percent of the true air ionization 
rate at 3 ft above a plane source of fission products uniformly distributed on 
the area. Therefore to obtain the true air ionization, in calculating exposure 
doses, the I values should be multiplied by 1.33. 

The average value of the terrain shielding factor, which is automatically 
contained in the source data, is about 0.75. The data, in general, apply to 
U-238 fission, for which the i fp (l) of Equation 3.65 is equal to 6.94 X 10~ 13 
R/hr at 1 hr per fission/sq ft or 3,510 R/hr at 1 hr per CT/sq mi. The value 
of i i <l) for the indicated induced activities is 0.13 x 10 13 R/hr at 1 hr 
per fission/sq ft. With these numbers, Equation 3.65 becomes 

k«T - 3.90 X 10- 13 [7717 * O.ow] «•“> 


K (1) 


1,970 r(l) + 0.019 


R/hr at 1 hr 
KT/sq mi 


(3.67) 


With Equation 3.67, the intensity-area integral becomes 

I g da = 1,970 [rU) + O.OI9] C(1)BW (3.68) 

It can be seen that, for unfractionated fission products (r(l) = 1) and 
for C(l) equal to unity, the ideal or upper limit value of the intensity-area 
integral per unit yield would be about 2,010 R/hr at 1 hr for observed fallout 
pattern intensities. 


31 






URS 702-1 


Values of the various fallout pattern features from the various scaling 
functions and Intensity-area Integral parameters are summarized In Table 3.4 
for selected weapon yields. The area integrals Indicate that the fraction of 
the total activity contained within the 1 R/'hr at 1 hr contour for the simpli¬ 
fied fallout pattern scaling systems varies from about 0.4 for a weapon yield 
of 1 KT to about 0.8 for a weapon yield of 100 WT. Although the scaling 
functions given In the previous sections apply only to weapon yields in the 
range of 30 KT to 100 NTT (for the latter high yield as an extrapolation), 
tabulated values for 1 and 10 KT are included in Table 3.4 as derived from 
the data of Reference 4. 

3.5 Ionization Rates and Potential Exposure Doses During Fallout Arrival 
3.5.1 Times of Fallout Arrival and Cessation 


The arrival time of fallout particles from stem altitudes Is 
estimated as the sum of two time periods. The first of these is the time 
required for the particles of a given average fall velocity to pass through 
the fireball and circulate in a toroidal motion until the updraft from the 
buoyant gases no longer carry the particles upward; the second time period 
is that of gravity fail under influence of winds to the ground. 

The apparent height at which particles of different diameters (and given 
fall velocities) are effectively free of the updraft under the rising fireball 
and cloud is estimated from 


Z = 1.80 X 10 V 


1 - exp (-0.00599W 


°- U6 .)] f<; 


. ^ 20 sec ; W =. 30 to 10 KT 


(3.69) 


The rate of increase of the apparent height with time after detonation at 
which the particles assume free fall is then given by 
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Table 3.4 

SUMMARY OF FALLOUT PATTERN FEATURES AND FALLOUT SCALING SYSTEM 
PARAMETER VALUES FOR AN ASSUMED EFFECTIVE WIND SPEED OF 15 MFH 


Pattern Feature 
or Quantity 

1 KT 

10 KT 

Weapon Yield 

100 KT 1 MT 

10 MT 

100 MT 

X 1 

-0.385 

-1.21 

-3.01 

-6.29 

-13.1 

-27.3 

*2 

0.0648 

0.342 

0.735 

0.882 

0.574 

-1.02 

*3 

0.352 

0.967 

2.06 

4.00 

7.68 

14.7 

X 4 

1.46 

5.18 

12.1 

19.3 

30.8 

48.6 

X 5 

0.834 

2,58 

5.00 

7.03 

4.39 

-14.3 

X 6 

1.34 

4.05 

8.56 

13.6 

21.6 

27,9 

X 7 

1.38 

5.30 

14.4 

28.8 

58.4 

118 

X 8 

1.91 

7.56 

20.8 

42.8 

88.6 

183 

x 9 

29.4 

61.2 

126 

258 

528 

1,080 

Y s (15) 

0.316 

0.796 

1.99 

5.02 

12.6 

31.6 

Y°(15) 

1.26 

2.28 

9.10 

33.6 

64.7 

123 

I 2,3 

1,730 

22,900 

35,300 

9,800 

3,550 

1,440 

*6 

130 

222 

540 

1,720 

5,270 

15,500 

I 7 

180 

317 

803 

2,712 

9,070 

30,510 

7 , b 







A /W 

93.4 

724 

544 

77.5 

14.4 

3.3 

a'/W 13 

757 

416 

687 

1,250 

1,412 

1,552 

A T /W b 

850 

1,140 

1,231 

1,328 

1,426 

1,555 

(7TT) + o ,oi9) c 

0.785 

0.787 

0.789 

0.793 

0.797 

0.803 

C(l) 

0.550 

0,735 

0.792 

0.850 

0.908 

0.983 

d 

z 2,3 

3,430 

8,910 

15,800 

22,900 

32,300 

44,400 

z d 

5,200 

15,600 

30,800 

40,200 

52,500 

68,500 

h d 

6,600 

18,400 

35,700 

52,300 

76,200 

111,000 

a Distances are in miles; 

intensit 

;ies are in 

R/hr at 1 

hr (observed); 

reference: 

100 percent 

fission 

yield 





b In R/hr at 1 hr per KT/sq ml; integrated to the 1 R/hr at 1 hr contour 
c r(l) calculated from data of Chapter 4 
d Altitudes in feet 
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Z f = 108 exp (-0.00599W~°‘ U6 t) ft/sec; t 2 20 sec; 

W = 30 to 10 5 KT (3.70) 

Approximate representations of the relationships between the particle fall 
rate, Z £ , at the height Z and the particle fall velocity vector, v £ , for a 
particle falling from the height Z to the ground (sea level) were derived for 
spherical particles from particle fall rate data as computed by Cassidy 15 ; 
these are; 

Z £ /v £ = 1.00 + 2.00 x 10~ 6 Z f ; v f < 2.0 ft/sec; Z = 0 to 10 5 ft (3.71) 

Z f /v f = 0.948 - 0.0530 log v £ + 3.78 X 10 _6 (1.259 log v £ )Z £ ; 

v f = 2.0 to 200 ft/sec; Z = 10 4 to 6 x 10 4 ft (3.72) 

The relatlonsnip among Z £ , v £ , and Z £ given by Equation 3.72 represents 
the tabulated data within 10 percent at the extremes of the indicated limits 
of application for v £ and Z. The indicated ranges in these parameters include 
all the values applicable to the fallout from the stem heights of interest ; 
therefore, Equation 3.72 generally is used, along with Equations 3.69 and 
3.70, to estimate the time after detonation when particles with a given fall 
vector leave the rising cloud and take on free fall. This time is given by 

t £ = 0.107W 0,116 |log[l585/v f + (1.259 + log v £ ) W 0-116 ] 

- log [o.779 (17.887 - log V £ ) + (1.259 + log v £ ) W° -116 ]|hr (3.73) 

The time for the particles with a given value of v £ to fall from the height 
Z £ (Z £ £ 60,000 ft) is estimated from 
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3.895 [2,034/v f + log v f - 17.887] 
[l,585w“ 


v f (1.259 + log v f) ] 

values between 2 and 
Equation 3.71, the ti 

t x = O.i.07*/ 0,116 [log (27.78 + W°* 116 ) - log (2,995/v f + W° ,U6 )]hr (3.7: 


Equations 3.73 and 3.74 apply to v f values between 2 and 200 ft/sec and to W 
values between 30 and 10^ KT. With Equation 3.71, the times are 


and 


9 [ 108 - 


(3.76) 


for v £ values less than 2 ft/sec and W values between 30 and 10° KT. 

The time of arrival on the ground, t & , for the particles with the falling 
velocity v f is equal to + t 2 * 

The particles that reach the earth's surface first are those with the 

falling velocity, v.- . For these particles the times after detonation 

2,3 

of leaving the cloud and of fall from Z f to the ground are given by 


5 [290.1W°* 0< 


(38.66 - log W)W 


f 0.116] 


log [o.77S (200.7 + log W) + (38.66 - log W)W°* 116 ]|hr; 


W = 30 to 10 KT 


and 


_i£SJLL 


Values of the earliest fallout arrival on the ground from stem heights as 
calculated from Equations 3.77 and 3.78, are given in Table 3.5 for several 
weapon yields. 
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The calculation of the time of fallout arrival and of fallout cessation 
(and the time period of deposition) for a given location on the stem fallout 
pattern canter line requires estimates of the minimum and maximum values of 
v f that arrive at the location. These two values of v f are for the particles 
that fall from the upwind and downwind sides of the stem volume; they are 
determined from the stem radius and various relationships between the down¬ 
wind distance, X, the fall vector, v^, and the apparent height, Z^, at which 
the particles take on free fall. 


Table 3.5 

ESTIMATED TIMES OF EARLIEST FALLOUT ARRIVAL FROM STEM ALTITUDES 


(hr) 

0.098 

0.167 

0.280 

0.463 


(hr) 

0.160 

0.260 

0.415 

0.655 


The maximum and minimum downwind distances of deposition along the center 
line of the stem fallout pattern for the particles that fall freely from the 
height, Z f , are estimated from 

X = 2.78 X 10 4 (v^/v f )Z f ± 0.0395W 0,333 exp (l.48 x lo" 4 w"°' 131 Z f ) miles; 

W ^ 30 to 10 5 KT (3.79) 

where v^ is in mi/hr; v^ is in ft/sec; and Z f is in feet. The value of Z f for 
a selected value of v f is estimated from 

10 6 (108 - v f ) 

1 ( 2 - 0v f + 5.990W -0 ' 116 


^ ft; v f < 2 ft/sec; W = 30 to 10 KT (3.80) 
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1.40 X 10 4 [2,034 - (17.887 - log v f )v f ] 
[(1.259 + log v f )v ( + 1.585W -0,116 ] 


ft ; 


v { = 2 to 200 ft/sec; W = 30 to 10 5 KT (3.81) 

The minimum and maximum values of v f are determined by selecting a series 
of values of starting with v f (2,3) as the largest value and then computing 
Z f and the two values of X in Equation 3.79; the v f values ara plotted as a 
function of X to determine the paired values of v f at a given value of X. 

Such a plot is shown in Figure 3.3 for the fallout from stem heights of a 5-WT 
yield detonation. The time of arrival of the particles with values from 
the minimum to the maximum, for a given value of X, is computed by substituting 
selected values of v f in Equations 3.73 and 3,74 or Equations 3.75 and 3,76 
and solving for both t^ and . 

Estimates of the v f values and arrival times of the particles deposited 
at locations which are a distance of V miles from the center line of the stem 
fallout pattern can be made by replacing the second term of Equation 3.79 
with 

± [l.56 X 10’V- 667 exp (2.96 X 10“V C,131 Z f ) - Y 2 ] / (3.82) 

The particles falling from cloud altitudes, the time of arrival on the 
ground of particles with the falling velocity, v f , may be generally represented 
by 


t B = 0.17W~° + t f hr (3.83) 

-0 i 

where 0.17W hr is the apparent time to cloud stabilization and t f is the 

time of fall from the height Z. The application of Equation 3.83 is discussed 
below in relation to: (1) the estimation of the time of arrival and cessation 
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of fallout at a given coordinate point, X, Y, with an assumed reference pattern 
spread due to lateral wind shear; (2) the estimation of the range in falling 
speeds (and diameters) of particles that arrive at a given location (as needed 
for estimating the gross activity-particle size distribution); and (3) the 
estimation of the influence of lateral wind shear ut the cloud heights on the 
spread of the fallout pattern. In all treatments, it is assumed that the cloud 
radius, r, remains constant in the direction of the wind (i.e., the continuing 
lateral expansion of the cloud after 0.17W 0-1 hr is neglected). 


The earliest arrival time of particles falling from the so-called stabi¬ 
lized cloud is that of the largest particles, with a fall vector designated 
by v^, from the lower surface of the cloud. The general area on which these 
particles fall is defined by an ellipse whose axis on the center line of the 
pattern is equal to r and whose axis in the cross-wind direction is equal to 
a distance designated as V . The first arrival time at downwind locations 
outside this elliptically shaped area of particles with falling speeds less 
than v f is for the particles falling from the forward edge of the cloud (at 
the cloud center height, h). The general equation for t f of the particles 
falling from the lower surface of the cloud that land at the point X, Y, is 


vO- -V(v— 


hr^b" 2 = 9.75W 0 ' 426 mi; W * 30 to 10 5 KT 

r = 0.464W 0 ' 431 mi; W = 1 to 10 5 KT 

r 2 b~ 2 = 3.07W -0 ' 262 mi; W = 1 to 10^ KT 

h” 1 = 0.314W -0,164 mi" 1 ; W = 30 to 10 5 KT 

2 -2 -O 272 w 

h b = 144W ■ ; W = 30 to 10' KT 


for 


and v^ in mi/hr. 
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The assumed equation for Y , the peripheral limit of th.' fallout pattern 
in the cross-wind direction is 

Y q = 2Y g (v w ) £l - e' k c (X c + r) ]miles (3.85) 

where Y g (v^) is given by Equations 3.41 through 3.44; this form of represen¬ 
tation assumes persistence of lateral shear c-u-nh that the pattern spread at 
very large downwind distances is twice the width, Y (v^), at the downwind 
distance, X g . Further, it is assumed that Y q is equal to 1.1 Y g (v w ) at x g 
so that 

= 0.8/(X b + r) miles -1 (3.86) 


« ' (,. . 1.395»°' 116 ) 


Substitution of v from Table 3.1 for v in Equation 3.84 givi 


(2. 9 , 0 - 0 < • v.x) 

“ 1,550W°'X* 8 , ' 

0.464W 0 * 431 ^(i.ssow 0 ^ 18 * v^)(i- y 2 /v^- 144(7.08 X W~°- 186 -vJ) 2 W 
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£l67W 

896W 


0.336 


0.082 



0.464W 0 ' 431 V(896W°-° 82 , V ;,)(l - vV)- 144(5.88 X vT 0 ' 254 - 


; W a 30 to 10 KT 


The approximate limits of X on the pattern center line (Y = 0) for Equation 
3.88 is X R ± r, and the approximate limits of X in Equation 3.89 is X g + r ± r. 
The hortest arrival time is for particles from the lowest point (center) of 
the cloud bottom, given by (h - b)/v g , which is 


or 


0.254 


: 0,186w“ , ‘“” (, - 0.0833W 0 - 136 ) hr; v < 3.28W u -'‘‘*° (3.91) 

The location on the pattern center line of earliest arrival is given by 


„.0.246 


1 f 


(3.92) 


The earliest particle fall time at downwind locations outside the area of 
the limiting ellipses defined by Equations 3.98 and 3.89 is "iven by 


t f - x c /v w 
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where X c is the downwind distance along the pattern center line to the center 
of a constant arrival time ellipse. The distance to the center of the (half) 
ellipse for an arbitrary point X, Y in the fallout area is defined by 

X = X - r%l - Y 2 /Y 2 miles (3.94) 

c ’ o 

In general, X must be solved by an iterative process using Equations 3.94 and 
3.85. The number of iterations required to obtain a solution of X c and Y q to 
within a few percent of the converged values is reduced if the first approxi¬ 
mation of in Equation 3.85. is taken to be equal to X - r if Y is less than 

r, equal to X if Y is more than 2r, and equal to X - 0.5r if Y is between r 
and 2r. If the Iteration is done by computer, the last given first approxima¬ 
tion of X could be used for all Y values, 
c 

The time for cessation (or last time of arrival) for all locations in 
the cloud fallout pattern is estimated from Equations 3.83 and 3.93 with X 
estimated from 

X c = X + r^l - Y 2 /Y 2 miles (3.95) 

together with Equation 3.85. If the center point of the constant arrival time 

ellipse is selected, the value of Y from Equation '*.85 is obtained directly 
o 

and the arrival and cessation time contours are then calculated from Equations 
3.94 and 3.95. 

3.5.2 L_ oendenc* of Particle Fall Times and Other Parameters on 
La .eral W1 \d Shear 


The particle fall times and scaling system parameters for the 
larger particles th»t fall to the ground from stem heights are assumed to be 
unaffected by lateral wind shear; thus lateral wind shear effects are considered 
only for particles falling from cloud heights. 
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The methods described below lor estimating the particle falling speeds 
and times and other parameters are for very simple lateral wind shear patterns. 
The basic assumptions for developing the mathematical representations are: 

(1) the lateral wind shear effect occurs only through the height range of the 
so-called stabilized cloud, i.e., from (h - b) to (h + b) ; (2) below the 
height, (h - b), all particles fall through the same set of meteorological 
conditions with a given average wind speed in the x direction, (thus the fall 
vectors of all particle groups for heights less than the bottom of the cloud, 
are parallel to the x axis); (3) the lateral wind shear through the height 
range of the cloud is either described by a given angle, 6 , to the average 
direction of the wind below the cloud height or by a single valued parameter, 
designated as cp, in degrees per unit increase in height above the bottom of 
the cloud; (4) the magnitude of the fall vector for each group of particles 
remains constant from the height of origin in the cloud to the ground; and 
(5) the magnitude of the wind speed is constant at all altitudes (a minor 
assumption to simplify the mathematical representations). 

The center of the cloud is taken as the center of the coordinate system; 
the periphery of the cloud is defined then by 


y ) - r (b - z*> = 0 


(3.96) 


The time of fall of particles with a given value of v £ from a point x, 
v, z j.?i the cloud to the bottom of tne cloud is given by 

t b (z + b)/v £ (3.97) 




K 


(3.98) 


where X b> Y b is the location of the particles in the z plane at (h - b) at the 
time, t . The lateral displacements are thus given by 
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X b - x = (v^/v f )(z + b)cos9 


(3.99) 


and 


Y b - y = (v w /v f )(z + b)sin9 


(3.100) 




6 = (z + b)cp degrees 


(3.101) 


or. If 9 Is constant, 


9=9 degrees 


(3.102) 


The time for the particle group to fall from the bottom of the cloud to 
the ground is given by 


t g = (h - b)/v f (3.103) 

or 

t = (X - XJ/v (3.104) 

g b w 

where X, Y, (Y = Y.) is the location at which the particles arrived on the 
b 

ground and X is given by 


x = x b + <V v f )(h " b) 


or, in combination with Equation 3.99, by 


X = x + (v w /v^)[(h - b) + (z + b)cos9] 


(3.105) 


(3.106) 
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t = (z + h)/v 


t = rex - x r + (y - y>i 


(3.107) 


(3.108) 


It should be noted that combinations of the above equations do not provide 
explicit solutions for cloud periphery values of x, y, and of z in terras of 
v^ and the cloud geometry equation constants when S depends on z as given by 
Equation 3.101. A solution for v f , however, is possible. (It may be noted 
that particles with the same value of originating at different heights in 
the cloud do not have the same diameter; for particles with a given diameter, 
v f is also a function of z.) 

Estimates of the fall velocity of particles originating from the top and 
bottom of the cloud volume and landing at the point X, Y may be made from 


[(h - b)X + (z + b)(Xcos9 + YsinO) ±Wr (1 - z /b K 


ill! 


where 


- [(h - b) + (z + b)cos9] 2 + (z + b) 2 sin 2 0 


(3.110) 


and 


H 2 = (h - b)Y - (z + b)(XsinS - YcosS) (3.111) 

for z values between -b and b. The limiting, or boundary, values among X, Y, 
and z at which real (as opposed to Imaginary) values of v f exist are obtained 
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by setting the quantities within the radical of Equation 3.109 equal to zero. 
The deposition pattern boundary lines for particle groups originating at the 
height z + h at the periphery of the cloud ellipsoid are thus given by 


<z + b)Xsin9 ± r^(l - z 2 /^ 2 )^ 
(h - b) + (z + b)cos9 


(3.112) 


The corresponding values of the fall vector mny be estimated fr< 
|[(h - b) + (z + b)cos9] X + (z + b)Y sin9 ) 


,(m) = - 


■ Y 2 - r 2 (l - z 2 /b 2 ) 


For a given value of X, Equation 3.112 represents a skewed elliptically- 
shaped curve in the Y , z coordinate system. The absolute maximum and minimum 
values of Y^ depend o.n both X and z. Although the slope of the lines in the 
X, Y plane increases as z increases from -b to b, the value of the intercept 
at X equal to zero has a maximum for particles that fall from the center height 
of the cloud (z = 0). Thus, if no restriction on v^(m) is given and 9 is posi¬ 
tive, the largest variation m would be 2r at X equal to zero and would 
approach the value 2bXsin9/[(h - b) + 2bcos9] at very large values of X. 

The general effect of lateral wind shear on the boundaries of the fallout 
patterns and on arrival and cessation times can be investigated through 
Equations 3.109 through 3.113. Although 9 could have any value between zero 
and ±180 degrees, the average lateral wind shear at heights from 40,000 to 80,000 
feet has been shown to be correlated with wind speed 12 whereby approximate values 
of Cp can be estimated from 

Cp = 108 v^ 1,14 degrees/mile (3.114) 

for the summer season, and 

cp = 38.3 v degrees/mile (3.115) 
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for the winter season; the applicable range of the two equations Is for z + h 
values In the range of about 7 to 15 miles. Thus, for an Integrated wind speed 
of 20 mi/hr, the average lateral shear angle from the bottom to the top of the 
cloud for a 1 WT yield explosion with a cloud thickness of 4.2 Miles would be 
about 15 for the summer season and about 30 for the winter season. The data 
on which Equations 3.114 and 3.115 are based are limited and perhaps somewhat 
outdated; further studies of similar data would be desirable to determine the 
frequency of occurrence of lar;er lateral wind shear angles than the average 
angles obtained from the above equations. 

If 9 is assumed or taken to be independent of z, the following solutions 
of z, x, and y, are obtained (the solution for is the same as that given by 
Equation 3.109 with 9 replaced by 9 ). 


v w b 2 j(Ysin9 o + X c os9 r ) v f - v w [t> + (h- b)cos9 o j| 
(r 2 v?*vV) “ 



(3.116) 
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where 


C 3 = r2y f + V w b [ (Ysln9 o + Xcos9 0 )v f ' v w (h ‘ b)cos9 o ] (3.119) 


and 


, 4 - [v,x - V „ - »]* 


(3.120) 


The limiting lateral extent of the fallout deposition, from the above 
equations, is defined by the condition that the Az, or Ax and Ay segments 
within the cloud be equal to zero (which is the same as the limit for obtaining 
real values of x, y, and z; the limiting value of Y for a given value of v f 
and X (with stated values of v and 9 ) is given by 


v bsin9 |r 2 v 2 + v bcos9 Txv - v (h - b)~]} 

Y = -2_ 21 _ t _ w oL f w_'Jl 

■ jxv, - ^- >(, - t 


An explicit equation for the limiting values of v f for the fallout at a 
given location is not possible since the limiting equation contains terms of 
v f with exponents from zero to 4. The values of v f that give maximum and 
minimum values of Y^ at any value of X may be estimated from 


/“(max) = (v /XH(h - b) + 


r 2 (h - b)b ■ + bcos9 ) 2 

+ —---- - -; 

2r (h - b)b cos9 + (X + bcos9 )‘ 
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(mlle3) (miles) 

10 10.3 

20 11.5 

50 15.5 

100 23.2 

500 94.0 


min 

niles) 

- 8.1 

-7 .2 

-4.9 

-3.0 

- 0.1 


AY° 

m 

(miles) 

18.4 
18.7 

20.4 
26.2 
94.1 


AY m /2r 

1.01 

1.03 

1.12 

1.44 

5.17 


The overall pattern spread at 100 miles away from ground zero for 30 of shear 
results in only a 44 percent Increase in the apparent pattern width relative to 
the "stabilized" cloud diameter. Although the center line of the pattern is 
slightly curved, the angle between the center line and the X axis to a distance 
of 50 or 100 miles for the illustration, is about 6°. 


3.5.3 Dependence of Particle Fall Times and Their Parameters on 
Total Wind Shear 


For overall wind shear effects, variation in wind speed as well as 
direction over the range of cloud thickness must be taken into account. To 
illustrate, the very simple shear condition is assumed in which v is the inte¬ 
grated wind speed up to the height of the cloud bottom and in which the integrated 
wind speed increases at a constant rate of S miles/hr per mile at greater heights 
(a variation in speed with height may also occur at lower heights, which would 
be accounted for in the integrated speed to the cloud bottom; this integrated 
wind speed would be approximately the same for all v^ if time variations in 
the wind speeds and directions at all z are not considered). The angle between 
the direction of the integrated wind speed vector to the bottom of the cloud 
and the direction of the wind above the bottom of the cloud is designated as 6. 

The direction of the integrated wind speed vector to the bottom of the cloud is 
defined to be parallel to the X axis. With these definitions and assumptions, 
the wind speed over the height range of the cloud is represented by 

v^(z) = v w + (z + b)S (3.127) 
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With Equation 3.127, an explicit solution for x, y, and z cannot be 
obtained for the Riven cloud geometry; however, the values of v^ for particle 
groups originating at the height z + h on the cloud periphery and landing at 
the point X, Y may be estimated from 


' V 


b)s] (z + b) (XcosQ - 


’ * * /»» >** - 


(3.128) 


wliere 


y 5 = | v w <h - b) + [v^ + (z + b)sj(z b)cos0j 
+ + (z + b)sj 2 (z + b)“sin 2 0 


(3.129) 


and 


V 6 = v w (h - b)Y - [v w + (z + b)sj (z + b) (XsinG - Ycos9) (3.130) 


The maximum values of Y for the particles originating at the height (z + h) and 
landing at the downwind distance X are given by 


|v w + (z + b)Sj(z + b)XsinS ± r V(’ - V 2 )-S 
v # (li - b) + [v + (z -r b)s](z + b)cos0 


(3.131) 


Although an explicit equation for the maximum and minimum values of Y is 
nut readily obtained, approximate values can be derived graphically. The effect 
of total shear on the spread of the pattern is illustrated by the following 






URS 702-1 


Y° calculated for a 

and an S value of 2. 

1 WT yield c 

. 1 per hour: 

loud, a 9 \ 

•alue of 30 

X 

(miles) 

Y° 

max 

(miles) 

'min 

(miles) 

£Y° 

(miles) 

<iY°/2r 

10 

10.6 

-7.9 

18.5 

1.02 

20 

12.3 

-6.8 

19.1 

1.05 

50 

17.8 

-4.6 

22.4 

1.23 

100 

28.4 

-2.8 

31.2 

1.71 

500 

124.7 

-0.1 

124.8 

6.85 


The overall pattern spread at 100 miles away from ground zero for 30 of 
lateral shear and an S value of 2.1 per hour of height above the bottom of the 
cloud results in a 171 percent Increase in the apparent pattern width relative 
to the "stabilized" cloud diameter. The average angle between the center line 
of the pattern and the X axis at a downwind distance of 100 miles is slightly 
larger than 7°. The values of for the particle group that land at X, Y m 
are given by the first term of Equation 3.128. 

The lateral displacement of particles with a given value of v^ originating 
at the height z + h on the cloud periphery and landing at the downwind distance 
X may be illustrated from 


Y 


(v w + (z + b)sj(z + b)sin9 




- b)s] 


The limiting value of v f for all real values of Y in terms of z and X are 
given by 
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v (li ~ b) + J"v *r C z r b)S*](z 1 b)cos6 
v (X) - —-^- - 1 - (3.133) 

X ± r VT"- z 2 /b 2 

nd the corresponding lateral displacement limit' i rom Equal ion 3.132, are 

fv + <7 + %J*| <r * bjsin'l , X ± r^l - z 2 /b s i 

y<x > • L Vii - 4 T«---b>8]<« ? r> o.j H «•***> 


The equations for v (X) and Y(X) give the envelopes of these two variables 
at a given value of X for all particle groups within the cloud volume. The 
paired values of each are for the particles originating from the cloud ellipse 
at y equal to zero. The general shapes of the v- envelopes are shown in 
Figure 3.4 for X equal to 50 and 100 miles. In addition, specific variations 
of with z + b are shown for the X, Y point.-, of 50,10 and 100,15. For each 
nf these curves, the values of z + b and v f that give the minimum and maximum 
of v f and t f are indicated; the points of tangenev with the envelopes are also 
shown as well as the major parameter values used ir. the calculations. 

An explicit general solution for z in terms of other parameters is not 
possible; the general equation of z is 

b 2 sV‘ + 2bS(v^b + 2S)z 3 

J- t jj - 2b"”S (XcosO + Ysiitf)j+ v w b-|v w + 2b v .. u)cosbj + Sb^Sfv^ + bS )J z“ 

* 2b 2 (v w + 2bS)£b(v w + bS) + v^(h - b)cos9 - v f (Xcos9 -t Ysin9)Jz 
+ b 3 (v^ + bS)|b(v w + bS) - 2^v f (Xcos0 + YsinS) - v^(h - b)cos0jj 

+ b 2 [xv f - v w (h - b)] 2 + v 2 fY 2 - r 2 ) = 0 (3.135) 

Approximate representations of the fallout pattern boundaries and the 
general directions of the boundaries, range of v values, and times of arrival 
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bo given in relation to the particle tra- 
c cloud originating at a given height along 
ular disc enciosj .g the area covered by the 


cles. The latter in given by 


(X - X) + (Y - V ) = r (1 - z /b ) 


+ (z + b)SJ (2 + b)cos5 


[v w + (z t b)sj(z + b)sin9 


The paired values of X^ and Y^, or the line on the X, Y plane along which the 
particles with any value (originating at the height z + b at the center of the 
cloud) are deposited, is given by 

[v w + (z + b)sj(z + b)X^sin9 

Y c = v w (h - b) + [v^ T (z * b)S](z + b)cosS (3.139) 

The values of z for particles with a given range in v f values landing at 
the location X , Y , from rearrangement of Equation 3.139 may be estimated from 


In Equations 3.137 to 3.139, the angle 9 may be a function of z; however, the 
separation of the two variables in Equation 3.140 indicates the assumption 
that 9 is independent of z. The value of v for the particle disc centered 


. X , Y , and origination at the height 


1 the cloud is represented 1 
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v^(h - b)sin9 
= X sin? - y _ cos9 


The value of Y , from combination of the above equations with Equation 
3.138, in terms of X and Y is given by 


r2 i 

Y = Y i (r/X )%X + Y 

me c 1 c 


The time of fall from the height z + h of the particle dis< 
X , Y , may be estimated from Equation 3,107 or from 


(X sin9 - Y sin9 'l 
V c o c ot 
2v S(h - b)sin5 


J 7 ~ - v . <h - b>y c 

f \ /x^slnti^ - Y c coj 


It may be noted that while Equation 3.136 defines the area covered by a 
circular disc in the X, Y plane centered at the location X^, Y , it defines a 
skewed elliptically shaped curve in the X fi , Y^ plane having a pseudo-center at 
the location X, Y; such a curve represents the loci of the particle disc 
centers on the X^, Y p plane for those discs whose radius is just equal to the 
distance between the disc center and the location X, Y. Thus the particles 

with the largest and smallest values of v ^ and t^ that deposit at the location 

X, Y are those on the edge of four different particle discs that originate at 
four different heights and whose centers are at four different X^ , Y^ locations. 
The values of X for particles landing at X, Y originating at the height z + h 
(and located at the edge of the particle disc) is given by 

Jh - b) + [v^ + (z + b)SJ(z + b)cosQ| 


±l: 


b)s"|(z + b)(Xcos9 + YsinS) 


■ b)s](z 
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The corresponding values of Y may be calculated from Equation 3.139. 

The variation of Y with X for the cited conditions for the X, Y locations 

50,10 and 100,15 is shown in Figure 3.5. The X , Y locations of the particle 

discs that are first and last to arrive (designated t g and t , respectively), 

that originate at the cloud mid-height (z = 0), and that originate from the 

maximum and minimum heights in the cloud (7 and Z , respectively) are 

max min 

indicated on each of the two curves. The X c , Y^ locations of the discs giving 
the maximum and minimum values of v f and are obtained graphically from a 
plot of t t as a function of v f . The respective values of these parameters are 
as follows: 


X 

(miles) 

50.39 

'55.99 

58.51 

42.53 

44.91 

49.49 


X = 50. Y 


10 


2.74 

4.37 

6.96 

8.50 

10.32 

12.09 


2.552 

2.406 

2.508 

4.191 

4.309 

4.081 


3.37 

3.73 

3.83 

2.63 

2.70 

2.92 


(z + b) 
(miles) 

0.836 

1.200 

1.850 
3.262 

3.850 
4,154 


100.56 7.07 

104.98 8.39 

108.55 12.25 

92.04 17.13 

95.23 20.06 

98.62 21.36 


X = 100, Y = 15 

1.319 
1.288 
1.334 
1.868 
1 .926 
1.890 


6.71 1.080 
6.98 1.225 
7.13 1.750 
5.77 3.003 
5.83 3.468 
6.00 3.584 


No explicit equations have been derived for estimating the maximum and 
minimum values of v f , t f , and z for the particles deposited at an arbitrarily 
selected X, Y location. 
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3.5.-1 Wind Speed and Shear Definitions for the Fallout Pattern 
Scaling Systems 

4 

As previously discussed, the fallout pattern scaling system 
coefficients were, to a large degree, evaluated from data on the fallout 
distributions from test shots Jangle S and Castle Bravo. Relative wind data 
for these two detonations arc applied here to Illustrate some definitions of 
wind speed and shear that are applicable to the scaling system. 

In general, the integrated wind speed refers to the total displacement 
of an imaginary particle of undefined size falling with constant speed from 
a given height (or altitude) to another height or to the ground. Thus, an 
integrated wind speed may be directly obtained from balloon sounding data. 
Where the latter are given as a function of azimuth and speed for given 
altitude increments, then v°(Z), the integrated wind speed from the ground 
surface to the height, Z, is calculated from 


2 w iri, i 


(3.145) 


where £Z^ is the thickness of the ith altitude increment, is the wind speed 
of the itn altitude increment, Z is the height or altitude of the nth increment, 
and 9° is obtained from 



2 V° s9 i 


i = 1 


(3.146) 


The average (unidirectional) wind speed is, in similar units, defined by 
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Tlit? integrated (or avvrage) particle displacement speeds for particles 
falling from large heights are usually less than the integrated wind speeds 
since solid particles (and especially those with large diameters) fali faster 
at the higher altitudes than at the lower altitudes. The integrated particle 
displacement speeds arc given by 


X v i Ai i sin9 i 



(3.148) 


where Lt ^ is the time of fall through the ith altitude increment, is the 
total time of fall through the n altitude increments, and is obtained from 



i = 1 


(3.149) 


both £ 1 t i and t^, for a given altitude increment and total distance of fall, 

\ary with particle diameter. 

If the wind speeds and directions remain constant over the period of fall 
of all the particles of interest from all Z, then a single value of v (Z) and 
v° for each Z is obtained. This condition may be approximately met for clouds 
from small yield detonations and for t n values of less than 2 to 3 hours. Eut 
for clouds from large yield detonations where t is greater than 2 to 3 hours 
even for fairly large particles, it is much less likely that the wind speeds 
and directions will remain constant in time at even one location; they generally 
are also not likely to be constant over the deposition area because of the 
resulting larger displacement distances. 
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The results from displacement calculations (details to be given in a 
following report) for the Jangle S Shot and Castle Bravo Shot clouds are 
summarized In Tables 3.0 and 3.7, respectively. 

The computed values ol v (/.), v (7.) for particle diameters from 75 to 
■100 microns and v' arc all about the same for the Jangle S Shot cloud; also, 
the 5 values differ very little from 6° for the integrated wind speed. 

However, because of the higher altitudes, variation of the wind sp eds and 
azimuths with lime (and space), and large changes in angle with altitude, 
the values of v^(Z), v^(Z), and for the Castle Bravo Shot show much larger 
differences; the same is true for 9^ and 9 n . The maximum difference in 5^ is 
8.5 for the particles from the Jangle S cloud and 2. < for the 1 0-micron 
particles from the Castle Bravo cloud. 

The accumulated (or integrated) shear factors relative to the height of 
the bottom of the cloud and the angular displacements for particles of different 
diameters are summarized in Tables 3.8 and 3.9 for the Jangle S and Castle Bravo 
8hot clouds, respectively. The total shear factor for the Jangle S Shot cloud 
is quite large, mainly due to the increase in wind speed with height above the 
cloud bottom; in addition, the values of both S and 6 are about the same for 
all the particle diameters from 75 to 400 microns. The values of S and 9 for 
the particles with 100- and 200-micron diameters from the Castle Bravo cloud 
include time and space variation; the variation in 9 with height is the major 
source of shear since the displacement of the particles starting at the top of 
the cloud is almost 180 from the directions of the hodograph of particles 
falling from the bottom of the cloud to the ground. 

The above described angular and total shear effects are included in the 
fallout pattern scaling functions. This fact should be kept in mind in applying 
the functions; in no case are the scaling functions directly applicable to a 
no-wind shear condition (t.e., the patterns always will contain residual effects 
of this shear, irrespective of the values of v_ used in the calculations). If 
a single wind speed and direction is assumed for all altitudes (a condition 
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Table 3.6 

SUMMARY OF INTEGRATED WIND AND PARTICLE DISPLACEMENT 
SPEEDS AND DIRECTIONS FOR THE JANGLE S SHOT CLOUD 


Wind or Particle Z - h , <1G 3 ft) 

Displacement - 

0 0,11 _ 1.11 _ 1 .48 2,11 2,96 

v°(7.) or v^(Z), (mi/hr) 


25.1 25.4 27.8 28.7 30.0 31.5 

24!8 25.1 27.4 28.3 29.5 30.9 


24.6 24.9 27.2 28.0 29.3 30.6 

24.6 24.9 27.2 28.0 29.2 30.6 

24.6 24.9 27.2 28.0 29.2 30.6 
24.5 24.8 27.1 27.9 29.1 30.5 

9° or 9^ (degrees) 

Integrated Wind 184.8 185.1 188.5 189.6 191.2 193.3 
Particle Displacement 

d = 75 microns 184.7 185.0 188.3 189.4 191.1 193.1 
d ^ 100 microns 184.7 185.0 188.2 189.4 191.0 193.1 
d = 200 microns 184.6 185.0 188.2 189.1 191.0 193.0 
d = 400 microns 184.6 184.9 188.2 189.3 190.9 193.0 


Average Wind Speed 
Integrated Wind Speed 
Particle Displacement 
Speed 

d = 75 microns 

d = 100 microns 
d = 200 microns 
d = 400 microns 




SUMMARY OF (RELATIVE) INTEGRATED r.'IND AND PARTICLE DISPLACEMENT 
SPEEDS AND DIRECTIONS FOR THE CASTLE BRAVO SHOT CLOUD 
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Table 3.8 

ACCUMULATED SHEAR FACTORS AND ANGLES FOR THE 
JANGLE S SHOT CLOUD 


Particle Diameter 
(microns) 


75 

100 

200 

400 


75 

100 

200 

400 


75 

100 

200 

400 


Z - h , (10 ft) 


v w (Z), (mi/hr) 


(24.6) a 

(24.6) 
(24.6) 
(24.5) 


41.8 

41.7 

41.9 

40.8 


42.5 

42.6 

42.7 
42.5 


43.3 

43.3 

43.4 
43.2 


43.7 

43.8 
43.8 
43.7 


6, (degrees) 


(184.7) a 15.0 

(184.7) 14.9 
(184.6) 14.8 
(184.6) 14.7 


15.9 16.6 
15.8 16.6 
15.8 16.6 
15.7 16.7 

(hr)" 1 


18.0 

18.0 

18.0 

18.0 


81.4 

81.5 
82.3 
8° .5 


63.8 

64.1 

64.4 

64.6 


46.6 

46.7 
46.9 
46.7 


34.1 

34.1 
34 .2 

34.2 


a Values in parenthesis are for the integrated particle displacements 
from h^. 
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Table 3.9 

ACCUMULATED SHEAR FACTORS AND ANCLES FOR THE 
CASTLE BRAVO CLOUD 


z - h b 
(IQ 3 ft) 

v w ( Z) 

(mi/hr 

) 

6 

(degrees) 

S 

(hr' 1 

)_ 


100 u 

200 w 

100 u 

200 u 

100 Li 

200 U 

0 

(19,6) a 

(18.2) a 

(254.8) 3 

(253.0) a 



3.75 

23.7 

17.9 

70.2 

30.1 

5.73 

-0.46 

C .75 

21.9 

16.3 

64.7 

69.6 

1.35 

-1.16 

13.75 

22.4 

16.9 

72.6 

78.0 

1.07 

-0.51 

18.75 

17.9 

13.1 

86.8 

96,6 

-0.48 

-1.44 

23.75 

15.7 

12.8 

114.3 

131.0 

-0.88 

-1.22 

25.06 

16.4 

13.6 

116.9 

132.9 

-0.68 

-0.98 

28.75 

18.3 

15.8 

124.3 

138.2 

-0.25 

-0.44 

33.75 

19.4 

17.6 

138.1 

151.7 

-0.031 

-0.10 

38.75 

22.3 

20.9 

148.5 

159.4 

0.37 

0.36 

43.75 

24 .0 

23.3 

168.6 

178.0 

0.53 

0.62 

48.75 

27.9 

27.0 

169.5 

177.8 

0.89 

0,95 

50.12 

29.2 

28.0 

170.8 

179.0 

1.00 

1.03 


Values in parenthesis are for the integrated particle displacements 
from h fa 
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that corresponds to the average unidirectional wind speed defined above), then 
the computed average wind speed must be reduced by a factor of 15/28.7 (or 
0.52) when using the scaling functions for yields less than 30 KT and by a 

factor of 15/26.6 (or 0.56) when using the scaling functions for yields larger 

* 

than 30 KT. These corrections are only approximate ones for the larger weapon 
yields; the correction factor actually should vary with particle size. However, 
these corrections will give a first order approximation of the pattern scaling 
function values when a single average wind speed is used. Revisions and exten¬ 
sions of the fallout pattern scaling system that are now in progress should 

minimize the wind speed correction problem; the major question to be resolved 
is whether the adjustments can be properly introduced without also introducing 
additional and more involved computational steps into the system. 

3.5.5 Distribution of the Condensed Radioactivity Among Particles as 
a Function of Falling Velocity Vector of Particle Diameter 

The activity carried by particles with a given value of v f falling 
from stem altitudes is conveniently represented by 

i^(l) = K f (l)A^ R/hr at 1 hr (3.150) 

where is the activity in fissions/ft 2 for a unit area of the stem (of speci¬ 

fied thickness) at the altitude from which the particles assume gravity fall, 

2 

and K f (l) is the conversion factor from flssions/ft to R/hr at 1 hr for the 
particle group. The total activity deposited at a given downvind distance in 
the stem fallout pattern is calculated from 

v f (max) 

I g = 1.467v w |(K f (l)A^/v 2 )dv f R/hr at 1 hr (3.151) 

v f (min) 


* No correction to the wind speed for the stem fallout pattern is suggested; for 
most applications a wind speed of 15 mi/hr should be used. 
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where v^(min) and v f (max) are the limiting fall velocity vectors in ft/scc 
for the smallest and largest particles that are deposited (as described in 
Section 3.5.1), and v is the average vxnd speed in mi/hr from the surface 
to Zj.. In the described scaling system for the stem fallout, only those 
particles with a given value of v arrive at :i given time; the arrival time 
for each particle group is given by the sum of Equations 3.73 and 3.74 (or 
3.75 and 3.76); thus the activity deposited up to a given time after detonation 
is calculated by graphical integration represented by 

Vj.(max) 

I g (t) = 3,38v^ | (K f (l)AVv*)dv f R. hr at 1 hr (3.152) 

Vf ( t) 

where v f (t) is the fall vector of the particle group arriving at t [v f (max) 
being the fall vector of the first particles that arrive]. The ionization 
rate at a given time is calculated by multiplying I g (t) b >' a decay correction 
factor, d f (t) so that 

I(t) = d f (t)I s <;) (3.153) 

where d f (t) is the average ratio of the ionization rate at the time, t, to the 
rate at H + 1 for all the particles that arc deposited at the location. If a 
large variation in the gross decay with particle sizes is to be assumed, then 
d f (t) could be applied in Equation 3.152 prior to the indicated integrations. 

To estimate the above integrals and the fallout deposition rates, the 
variation of ij(l) with v^ must first be computed from the stem fallout pattern 
scaling functions. This is done by using the variation of I with X, for 
X 2 (x^ + X^)/2 and v^ = 15 mph, along the center line of the stem fallout 
pattern and, in the first approximation, taking a single value of K f (l)A^ in 
Equation 3.151 for the particles that fall at a location. The single (average) 
value of K f (l)A^, is then written as 
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K f (l)A' = -=- - -=■ R/hr at 1 hr (3.154) 

1.467v Il/v f (min) - l/v f (max) 

The first estimate of the value of v^, for association with each K f (l)A^ value 
is computed from 

log v f = -0.050 + 0.769W 0,048 log v f (min); W =“ 10 3 to 5 x 10 4 KT (3.155) 


log v f = 0.070 + 1.25w“ 0,043 log v f (max); W “ 10 3 to 5 X 10 4 KT (3.156) 


or from the average of the v values computed from both equations. The adjusted 
2 1 

values of K f (l)A^/v f as a function of v f are determined from successive approx¬ 
imations of I g in Equation 3.151 starting with K f (1)Aj for K f (l)A^ and v f for v f . 

The activity deposited by particle groups with a given fall vector from 
cloud altitudes at a given downwind location, for the simple fallout pattern 
scaling system in which a uniform distribution of each particle throughout the 
cloud volume is assumed, is represented by 




(•:< 


T 


R/hr at 1 hr 


for locations on the fallout center line with either a directionally balanced 
effective (or no) lateral wind shear. In Equation 3.157, A is the average 

3 * 

initial cloud concentration in fisslons/ft of the activity carried by the 

particles with the fall vector v , K (1) is the conversion factor in R/hr at 
2 11 

1 hr per fission/ft , and the remaining combined terms represent the path 
length through the elliptical cloud from which the particles with the fall 
vector v^ originate and which deposit at the downwind distance X on the pattern 
center line. Integration of l f (l) over all the particle groups that land at 
the location gives the total activity deposited; this is represented by 
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v (max) 

l g = 1.467v^ I (i £ (l)/v £ )dv R/hr at 1 hr (3.158) 

v f (min) 

for In mi/hr and v £ In ft/scc. 

The average deposition rate of the particle group with the fall vector, 
v £ , from the cloud altitudes is given by 

i f (l)/&t. f = 3.6 X 10 3 v f K f (l)A f R/hr at 1 hr/hr (3.159) 

where At f is the time period over which particles with a given value of v land 
on the ground; £t £ may be obtained from a plot of t £ as a function of v f ; for 
locations near the pattern center line, for the case of no wind shear, it is 
given by 


4; - v-*) - h - S 

EW 


The deposition rate of all the particle groups arriving at a given time 
after detonation at a location is represented by 


(dl /dt) = 5.28 X 10‘ 


I 

.. | « f U)A f /, f )dv f 


R/hr at 1 hr/hr (3.161) 


Integration of Equation 3.161 gives, for I (t), 
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I (t) = 5.28 x 10 3 v ( dt 

s w 1 


The estimates of K f (l)A f for difference values of v f are made through a 

series of approximations starting with selected values of I and the downwind 

distance along the center line of the cloud fallout pattern at which they are 

calculated to occur, using the simplified scaling system with a 15 mph average 

wind speed. The following equations may be utilized to provide the first 

estimates of K f (l)A f ar.d v f (designated at K f (l)A f and v^, respectively) for 

each set of paired values of I and X. 

s 


f (max,t) 
K (1 )A 


- dv R/hr at 1 hr 


(3.163) 


where 


= (X + r) W r 2 h 2 b ~ 2 4 (X 7 r ) 2 

(X - r) +Vr 2 hV 2 + (X - r ) 2 


(3.164) 


and 


12L 


r) W r 2 h 2 b- 2 7 (X 7 r) 2 

x o+ VXv r 77 


(3.165) 


70 




URS 702-1 


in which 


r % 0.464W ’ miles, W = 1 to 10' KT 


X = 0.186v 1 


miles, W = 30 to 10 KT 


itimate of v for each value of K (1 )A . is obtained f r< 


: 30 to 10 KT 


for X in miles. If the estimates are made using wind speed other than 15 mi/hr, 
the coefficient 70.0 is replaced with 4.67v (for v in mi/hr). 

The potential exposure dose during the fallout arrival period is calculated 
by integrating graphically the computed values of I(t) from t^ to t for the 
fallout from both stem and cloud altitudes. Also, the value of K^(l) is 
adjusted to eliminate the instrument response factor, l, of equation 3,65; 
this correction, for the I g <t> values computed through use of the simplified 
fallout scaling system, requires increasing th< integrated doses by a factor 
of 1.33. 

Because the rather complicated form of the above equations giving the 
variation of I g (t) and I(t) with time after detonation, only a few example 
calculations were made from which approximating functions were derived to 
represent the variation of I g (t) or the accumulation of .he fallout from stem 
and cloud heights with time after detonation. Since the model of the stem 
depicts the fall of the larger particles from lower heights in combination 
with the gradual dilution of the particle concentrations with height (and 
time) and the deduced activity-size distribution (from the arrangement of the 
stem fallout pattern), a fallout rate given as dl g (t)/dt is obtained that is 
highest at the first arrival of fallout and decreases as an inverse function 
of time during the fallout period. On the other hand, the model for the 
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fallout from the cloud wherein the particles are assumed to fall from a 
stationary source volume of uniform concentration, the fallout rate increases 
from aero at the first arrivul of fallout, Increases to a maximum, and then 
decreases to zero again at fallout cessation. This general picture of the 
deposition rates is fairly well verified by analysis of observed intensity 
(exposure rate) data; however, detailed data from time-interval collectors 
show large short-time variabilities in the deposition rates which indicates 
that the particles are not uniformly distributed in the source volumes. 

(Some of the observed variability in the collector data is undoubtedly also 
caused by local variations in the surface wind speeds and directions over the 
deposition period and the apparent efficiency of the collectors.) 

The simplified versions of the increase of I g (t) with Vime during the 
fallout period are given in terms of the fraction of I g that is deposited at 
a given fraction of the deposit period; the representation is defined by 

I g (t) = f(t)I g (3.1G7) 

and the fraction of the time during the deposition period is defined by 

T = (t - t )/(t - t ) (3.168) 

a c a 

For the fallout from stem heights with f(t) redesignated as f_, the above- 
described variation of the fallout rates with time is represented approximately 
by 


t°Jl - (2r - l) 2 

df s /dT = ----- (3.169) 

where f g is a constant; integration of Equation 3.169 under the condition that 
f g is zero when T is zero and f g is equal to one at t equal to one gives 

i g = (1/tt)^ 1 - (2t - l) 2 + sln' l (2T - 1) j + 0.500 (3.170) 

where f is equal to 1 /tt. 
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For the fallout from cloud heights with f(t) redesignated as f , the 
above-described variation of the fallout rates with time is represented 
approximately by 


for a symmetrical rate of arrival curve that has a peak rate of urrival at the 
mid-time of fallout period (i.c. , at ' = 0.5), Integration oi Equation 3.171 
gives 



In regions of heavier fallout from the cloud heights, the peak fallout rate 
generally occurs somewhat earlier in the fallout period than at the mid-time, 
more nearly at the time of r equal to 0.4 than 0.5. The equations for f with 
the peak rate of fallout at - equal to 0,4 are 




The calculated variation of f and of f with t is given in Figures 3.6 

and 3,7, respectively. The curves can be used directly for estimating I(t); 

for the cloud fallout, and the values of f taken from the curves for T equal 

c o 

0.4 would be preferred. At the shorter downwind distances the peak fallout 
rate may occur even earlier in the fallout period; also, the time of the peak 
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ra^e would tend to shift towards the mid-time as the downwind distance 
becomes very large. The approximations for f g and f^ given by Equations 
3.170, 3.173, and 3.174 are consistent with the model formulation for 
application to computer programming in damage assessment evaluations. 

To illustrate the shapes of the I(t) curves as a function of time after 

detonation, values of this variable were calculated for the following assumed 

conditions; for the stem fallout, t is 0.5 hrs, t is 1.0 hr, I is 1,000 
a c s 

R/hr at 1 hr, and d(t) is equal to 1/t; for the cloud fallout, t is 1.0 hr, 

t is 2.0 hr, I is 1,000 R/hr at 1 hr, and d(t) is equal to 1/t. Since, in 
c s 

general, 

t = t + (t - t )t (3.175) 

a c a 

then 

t = 0.5 + 0.5T (3.176) 

for the assumed stem fallout case, and 

t = 1 + T (3.177) 

for the assumed cloud fallout case. Also, 

I(t) = 1,000 f g /t (3.178) 

for the stem fallout, and 

I(t) = 1,000 f /t (3.179) 

c 

for the cloud fallout. The calculated values of I(t) for the two assumed sets 
of conditions are plotted as a function of t in Figure 3.8. The peak intensity 
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for the assumed stem fallout case Is almost 1,100 R/hr and it occurs at about 
0.8 hrs (t = 0.6). The peak Intensity for the assumed cloud fallout case is 
505 R/hr and it occurs at about 1.9 hrs (t = 0.9). 

The air ionization rate from the airborne particles can be estimated 
from either the calculated deposition rates or the estimates of K f (l)A f . 

From calculations made by Laumetz,*'® the air ionization rate at three feet 
inside a plane boundary of a semi-infinite volume containing a uniformly- 
distributed radiation source is given by 

1 (t) = 0.950 X 10~Va (t)E, in R/hr (3.180) 

a tot 

in which A t is the activity in dis/sec per fission at the time, t, after 

fission, A (t) is the concentration in fissions/ft 3 , and E is the total 
o t 

gamma ray energy in Mev/dis. 

At a given time, A and E are constant, so that the contributions to 
activity in A f from all the particle groups arriving at a given time can be 
summed as in the case of the deposited material. Multiplication by the appro¬ 
priate values of A^ and E^ for the designated time will convert the sum to 
R/hr. The fallout from stem altitudes arriving at a given instant consists 
of particles with a single value of v^ whose air activity concentration is 
estimated from 

A o <0 = k f ( l jv f (t) F I ) fissions /ft 3 <3181) 

The values of A q for the fallout from cloud altitudes arriving at a given 
Instant is estimated from 

v f (max,t> 

A (t) = 3.38v 1 (A,/v )(d log v„) fissions/ft 3 (3.182) 

o w I f f f 

v f (min,t) 
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in which A f is equal to i f /K f (l) for the particles arriving at the time, 
and where is in mi/hr and v f is in ft/sec. 
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SUMMARY 


This report summarizes the fallout pattern scaling relationships that 
were developed in the period 1962 to 1964; the report Includes the values 
of the scaling equation coefficients that were derived from selected fallout 
pattern data. The meaning of the scalar wind speed multiplier that is used 
in the scaling equations is discussed relative to computer applications of 
the scaling system and approximate wind speed adjustment factors for use 
with wind speed averages that may be assumed in such applications are pro¬ 
vided. The relative degree of wind shear inherent in the scaling system 
parameters is also discussed in some detail. Basic equations for relating 
surface density of radionuclides and air ionization rates including consider¬ 
ation of fractionation, surface roughness, and instrument response are given 
and discussed together with the influence of these factors and others on the 
limiting values of K factors that represent the relative amount of the 
radioactive sources contained within the deduced area covered by the fallout 
patterns. Scaling equations and data are also presented for use in esti¬ 
mating, for any location in the fallout region, the time of fallout arrival, 
the time of fallout cessation, the variation of the exposure rate (i.e., air 
ionization rate in roentgens per hour) with time during fallout arrival, and 
the total potential exposure from the time of fallout arrival to selected 
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Appendix A 

SUMMARY OF SELECTED SCALING SYSTEM PARAMETERS 

Provisional values of several of the fallout scaling system model 
parameters for weapon yields In the rai’ge of 1 to 30 KT, not provided In 
the main text of the report, arc as follows: 

X 2 =" 0.0139W 0,194 - 10.34W"°' 167 J miles; W = 1 to 30 KT (A-l) 

X 3 “■ 0.0139W 0,494 [v w + 10.34W _0,157 J miles; W * 1 to 30 KT (A-2) 

X 4 **■ 0.0973v w W°' 550 miles; II = 1 to 30 KT (A-3) 

Z 2 3 * 3,43 X lo3w °' 414 feet : W = 1 to 30 KT (A-4) 

l 2 3 ~ 3,483 X W 1 ’ 122 : W = 1 to 30 KT (A-5) 

m w = 0.793; W = 1 to 30 KT (A-6) 

X x = -0.385W 0,496 miles; (v^ = 15 mi/hr); W = 1 to 30 KT (A-7) 

Xj * Xg - 0.139W 0,404 log Ig 3 ; W = 1 to 30 KT (A-8) 

X 5 = 0.0556v w W 0,46 ' miles; W = 1 to 30 KT (A-9) 

X g = 0.0894v w W° -481 miles; W = 1 to 30 KT CA-10) 

X ? =■ 0.0919v w W°' 586 miles; W = 1 to 30 KT (A-ll) 

X g = 0.128v w W° -596 miles; W = 1 to 30 KT (A-12) 

X g = 1.956v w W 0,319 miles; W = 1 to 30 KT (A-13) 

v^r/h = 8.34W -0 ' 03 mi/hr; W = 1 to 30 KT (A-14) 

“0.050 

v^r/h = 5.19W mi/hr; W = 1 to 30 KT (A-15) 

v ? r/h = 5.04W -0,155 mi/hr; W = 1 to 30 KT (A-16) 

A . 1 
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The values of K°, n^, (v 2 r 2 b 2 )°, and for the yield range of 1 to 30 KT 
are the same as those given in Table 3.1 for application in Equations 3.28 
and 3.30. 

Calculations of fallout arrival times for average wind speeds of 10, 20, 
and 30 mi/hr at selected downwind (X) and crosswind (Y) distances for weapon 
yields of 1, 3, and 10 MT are summarized ir. Tables 1A, 2A, and 3A, respectively. 
The wind speed corrections discussed in paragraph 3.5.4 do not apply to the 
arrival time calculations. Similar calculations on the fallout cessation times 
are summarized in Tables 4A, 5A, and 6A. 


Calculations of the exposure dose during the fallout period and from 
fallout cessation to H + 36 at selected X,Y locations are summarized in Table 
7A in terms of Exposure Dose-Standard Intensity Ratios for a detonation yield 
of 1 MT and an average wind speed of 20 mi/hr. The calculations show that 
D(36)/I g is essentially independent of Y at a given value of X, The values of 
D(36)/I g are represented very closely by 

D(36)/I g = 3.16 e" 9,79 X 10 X ; X S 85 km, All Y (A-17) 

and 

_3 

D(36)/I g = 1.73 e" 2,72 X 10 X ; X * 85 km, All Y (A-18) 


tor W = 1,000 KT, v^ = 20 mi/hr, and X in km. The additive exposure doses 
for times greater than H + 36 may be obtained from the so-called dose-rate 
multipliers (DRM) given in Volume II of Reference 4. The variation of I(t)/I 
with time after detonation for downwind distances of 10, 20, 50, 150, 250, and 
350 km is illustrated in Figures 1-A through 6-A. 


A.2 
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Table 1A 

SUMMARY OF TIME OF ARRIVAL OF FALLOUT FROM CLOUD HEIGHTS 
AT SELECTED X,Y LOCATIONS FOR W = 1,000 KT 
(Time in Hours After Detonation) 

± Y (km) 

X (km) 0 5 10 20 40 60 80 100 120 140 

1. v — 10 mi/hr 

10 0.94 0.94 0.94 0.96 1.12 

20 0.95 0.95 0.96 0.97 1.07 

30 1.04 1.04 1.06 1.12 1.36 

50 2.28 2.84 2.29 2.32 2.44 2.65 

150 8.50 8.50 8.50 8.51 8.55 8.61 8.71 8.87 9.14 

250 14.71 14.71 14.71 14.72 14,75 14.80 14.88 15.00 15.18 

350 20.93 20.93 20.93 20.94 20.96 21.01 21.09 21.19 21.36 21.71 

2. v =20 mi/hr 

10 0.49 0.49 0.51 

20 0.50 0.50 0.52 

30 0.56 0.58 0.68 0.81 

50 1.18 1.19 1.21 1.30 

150 4.29 4.29 4.30 4.32 4.40 4.60 

250 7.40 7.40 7.40 7.41 7.46 7.57 

350 10.51 10.51 10.51 10.52 10.56 10.64 10.81 

450 13.61 13.61 13.62 13.63 13.66 13.74 13.88 

550 16.72 16.72 16.72 16.73 16.77 16.84 16.97 

650 19,83 19.83 19,83 19.84 19.88 19.94 20.08 

750 22.94 22.94 22.93 22.95 22.98 23.05 23.18 

3. v =30 mi ''hr 


10 ----- - 
20 0.49 0.50 0.52 

30 0.51 0.52 0.55 

50 0.82 0.83 0.87 

150 2.89 2.89 2.90 2.93 3.09 

250 4.96 4.96 4.97 4.98 5.06 

350 7.02 7.03 7.04 7.05 7.10 7.24 

450 9.10 9.10 9.M 9.12 9.16 9.26 

550 11,18 11.18 11.13 11.19 11.23 11.32 

650 13.25 13.25 13.25 13,26 13.30 13.38 

750 15.32 15.32 15.32 15.33 15.37 15.45 

850 17.39 17.39 17.39 17.40 17.44 17.52 

950 19.46 19.46 19.46 19.47 19.51 19.59 

1050 21.53 21.54 21.54 21.54 21.58 21.66 

1150 23.61 23.61 23.61 23.62 23.65 23.73 


A .3 
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Table 2A 

SUMMARY OF TIME OF ARRIVAL OF FALLOUT FROM CLOUD HEIGHTS 
AT SELECTED X,Y LOCATIONS FOR W = 5,000 KT 

Y (km) 

X (km) 0 5 10 20 40 60 80 100 120 140 

1 . v w = 10 mi/hr 

10 1.16 1.17 1.17 1.19 1.27 

20 1.15 1.15 1.15 1.17 1.23 

30 1.20 1.20 1.20 1.21 1.28 

50 1.72 1.73 1.74 1.78 1.95 2.23 

150 7.94 7.94 7.94 7.95 7.99 8.07 8.18 8.33 8.54 8.86 

250 14.15 14.15 14.15 14.16 14.19 14.24 14.31 14.41 14.55 14.73 

350 20.37 20.37 20.37 20.38 20.40 20.44 20.51 20.59 20.71 20.86 

2. v =20 ml/hr 

10 

20 1.29 1.31 

30 1.15 1.16 1.17 1.22 

50 1.21 1.21 1.22 1.26 

150 4.01 4.01 4.01 4.04 4.14 4.36 

250 7.11 7.12 7,12 7.13 7.19 7.30 7.48 

350 10.22 10.22 10.23 10.24 10.28 10.36 10.48 10.72 

450 13.33 13.33 13.33 13.34 13.38 13.44 13.55 13.72 

550 16.44 16.44 16.44 16.45 16.48 16.54 16.63 16.79 

650 19.54 19.55 19.55 19.55 19.59 19.64 19.73 19.87 

750 22.65 22.65 22.65 22.66 22.69 22.75 22.83 22.97 23.27 

3. v =30 mi/hr 

10 0.58 0.59 0.65 

20 0.55 0.57 0.58 

30 0.56 0.56 0.58 

50 0.65 0.66 0.71 

150 2.70 2.70 2.71 2.75 2.95 

250 4.7'/ 4.77 4.77 4.79 4.88 5.09 

350 6.84 6.84 6.84 6.86 6.92 7.04 

450 8.91 8.91 8.92 8.93 8.97 9.06 9.26 

550 10.98 10.98 10.99 11.00 11.04 11.11 11.25 

650 13.06 13.06 13.06 13.07 13.10 13.17 13.29 

750 15.13 15.13 15.13 15.14 15.17 15.23 15.35 

850 17.20 17.20 17.20 17.21 17.24 17.30 17.40 

950 19.27 19.27 19.27 19,28 19.31 19.37 19.47 

1050 21.34 21.34 21.34 21.35 21.38 21.44 21.54 21.79 

1150 23.41 23.41 23.42 23.42 23.45 23.51 23.60 23.84 


A.4 




URS 702-1 


Table 3A 

SUMMARY OF TIME OF ARRIVAL OF FALLOUT FROM CLOUD HEIGHTS 
AT SELECTED X,Y LOCATIONS FOR W = 10,000 KT 

Y (km) 

X (km) 0 5 10 20 40 60 80 100 120 140 

1. = 10 mi/hr 

10 1.46 1.46 1.48 1.48 1.54 1.65 

20 1.44 1.44 1.44 1.45 1.50 1.59 

30 1.45 1.45 1.45 1.46 1.51 1.59 

50 1.57 1.57 1.58 1.59 1.64 1.75 

150 6.93 6.93 6.94 6.95 7.01 7.11 7.24 7.43 7.66 7.96 

250 13.15 13.15 13.15 13.16 13.19 13.25 13.33 13.43 13.57 13.74 

350 19.36 19.36 19.37 19.37 19.40 19.44 19,50 19.59 19.69 19.82 

2. = 20 ml/hr 

10 

20 1.64 1,65 1.69 

30 1.48 1.48 1.49 1.54 

50 1.44 1.44 1.45 1.48 1.65 

150 3.50 3.50 3.51 3.55 3.70 3.97 

250 6.61 6.61 6.61 6.63 6.71 6.84 7.04 7.39 

350 9.72 9.72 9.72 9.73 9.78 9.87 10.00 10.20 10.54 

450 12.82 12.82 12.83 12.84 12.68 12.94 13.05 13.20 13.42 

550 15.93 15.93 15.93 15.94 15.98 16.03 16.12 16.25 16.42 16.71 

650 19.04 19.04 19.04 19.05 19.08 19.13 19.21 19.32 19.48 19,71 

750 22.15 22.15 22.15 22.15 22.18 22.23 22.31 22.41 22.55 22,75 

3. v =30 mi/hr 


10 

20 

30 

50 1.51 1.52 1.54 

150 2.36 2.36 2.38 2.44 

250 4.43 4.43 4.44 4.47 

350 6.50 6.50 6.51 6.52 

450 8.57 8.57 8.58 8.59 

550 10.64 10.64 10.65 10.66 

650 12.71 12.72 12.72 12.73 

750 14.79 14.79 14.79 14.80 

850 16.86 16.86 16.86 16.87 

950 18.93 18.93 18.93 18.94 

1050 21.00 21.00 21.00 21.01 

1150 23.07 23.07 23.08 23.08 


2.71 

4.59 4.83 

6.60 6.74 7.00 

8.64 8.75 8.91 

10.70 10.78 10.91 11.14 

12.77 12.83 12.94 13.12 

14.83 14.89 14.99 15.14 15.43 

16.90 16.96 17.04 17.17 17.40 

18.97 19.02 19.10 19.22 19.42 

21.04 21.09 21.17 21.28 21.46 

23.11 23.16 23.23 23.34 23.51 
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Table 4A 

SUMMARY OF TIME OF CESSATION OF FALLOUT FROM CLOUD HEIGHTS 
AT SELECTED X,Y LOCATIONS FOR W = 1,000 KT 

Y (km) 

X Ocm) 0 5 10 20 40 60 CO 100 120 140 

1. v =10 mi/hr 

10 1.62 1.61 1,60 1.56 1.28 

20 2.24 2.24 2.23 2.20 2.05 

30 2.86 2.86 2.85 2.83 2.72 

50 4.10 4.10 4.10 4.08 4.01 3.87 

150 10.32 10.32 10.32 10.31 10.27 10.21 10.12 9.98 9.71 

250 16.53 16.53 16.53 16.52 16.50 16.44 16.37 16.25 16.08 

350 22.75 22.75 22.75 22.74 22.71 22.66 22.59 22.48 22.32 21.97 

2. v w = 20 mi/hr 

10 0.85 0.84 0.80 

20 1.16 1.15 1.13 

30 1.47 1.47 1.45 1.37 

50 2.09 2.09 2.08 2.03 

150 5.20 5.20 5.20 5.18 5.11 4.93 

250 8.31 8.31 8.31 8.29 8.25 8.15 

350 11.42 11,42 11.41 11.40 11.38 11.29 11.12 

450 14.52 14.52 14.52 14.51 14.48 14.40 14.26 

550 17.63 17.63 17.63 17.62 17.58 17.51 17.38 

650 20.74 20.74 20.74 20.73 20.69 20.62 20.49 

750 23.85 23.85 23.84 23.84 23.80 23.73 23.60 

3. = 30 mi/hr 

10 

20 0.80 0.79 0.74 

30 1.01 1.00 0.97 

50 1.42 1.42 1.40 

150 3.50 3.49 3.49 3.47 3.32 

250 5.57 5.57 5.56 5.55 5.48 

350 7.64 7.64 7.64 7.62 7.57 7.44 

450 9.71 9.71 9.71 9.70 9.65 9.56 

550 11.78 11.78 11.78 11.77 11.73 11.64 

650 13.85 13.85 13.85 13.84 13.80 13.72 

750 15.93 15.93 15.92 15.92 15.88 15.80 

850 18.00 18.00 18.00 17.99 17.95 17.87 

950 20.07 20.07 20.07 20.06 20.02 19.95 

1050 22.14 22,14 22.14 22.13 22.09 22.02 

1150 24.21 24.21 24.21 24.20 24.17 24.09 
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Table 5A 

SUMMARY OF TIME OF CESSATION OF FALLOUT FROM CLOUD HEIGHTS 
AT SELECTED X,Y LOCATIONS FOR W = 3,000 KT 

Y (km) 

X (km) 0 5 10 20 40 60 80 100 120 140 

1. v =10 mi/hr 

i.O 2.16 2.16 2.15 2,11 1.93 

20 2 78 2.78 2.77 2.74 2.61 

30 o.40 3.40 3.40 3.37 3.27 

50 4.65 4.65 4.54 4.62 4.55 4.41 

150 10.86 10.86 10.86 10.85 10.82 10.76 10.67 10.54 10.36 10.07 

250 17.08 '7.08 37.07 17.07 17.04 16.99 16.93 16.83 16.71 16.54 

350 23.29 23.29 23.29 23.28 23.26 23.22 23.15 23.07 22.96 22.81 

2. = 20 ml/hr 

10 1.11 

20 1.43 1.42 

30 1.74 1,73 1,72 1.C4 

50 2.36 2.36 2.35 2.30 

150 £.17 5.47 5.46 5.45 5.37 5.20 

250 8.58 8.58 8.57 8.56 8.51 8.42 6.25 

350 11.68 11.68 11.68 11.67 11,63 11.56 11.44 11.21 

450 14.79 14.79 14.79 14.78 14.75 14.68 14.58 14.41 

550 17.90 17.90 17.90 17.85 17.86 17.80 17.71 17.56 

650 21.01 21.01 21.00 21.00 20.97 20.91 20.82 20.68 

750 24,11 24.11 24.11 24.10 24.07 24.02 23.93 23.80 23.50 

3. = 30 ml/hr 

10 0.77 0.75 

20 0.98 0.97 

30 1.19 1.18 

50 1.60 1.59 

150 3.67 3.67 

250 5.74 5.74 

350 7.81 7.81 

450 9.89 9.89 

550 11.96 11.96 

650 14.03 14.03 

750 16.10 16.10 

850 18.17 18.17 

950 20.24 20.24 

1050 22.32 22.32 

1150 24.39 24.39 
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0.69 
0.92 
1.14 
1 .57 

3.66 3.63 3.49 

5.74 5.72 5.65 

7.81 7.80 7.75 

9.88 9.87 9.83 

11.96 11.95 11.91 

14.03 14.02 13.98 

16.10 16.09 16.06 

18.17 18.16 18.13 

20.24 20.23 20.20 

22.31 22.31 22.28 

24.39 24.38 24.35 


5.46 

7.64 

9.75 9.56 

11.84 11.70 

13.92 13.80 

16.00 15.89 

18.07 17.97 

20.15 20.05 

22.22 22.12 21.87 

24.29 24.20 23.97 
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Table 6A 

SUMMARY OF TIME OF CESSATION OF FALLOUT FROM CLOUD HEIGHTS 
AT SELECTED X,Y LOCATIONS FOR W = 10,000 KT 

Y (km) 

X (km) 0 5 10 20 40 60 80 100 120 140 

1. v w = 10 ml/hr 

10 3.15 3.14 3.14 3.11 2.98 2.71 

20 3.77 3.77 3.76 3.73 3.63 3.42 

30 4.39 4.39 4.38 4.36 4.27 4.10 

50 5.63 5.63 5.63 5.61 5.54 5.41 

150 11.85 11.85 11.84 11.84 11.80 11.74 11.65 11.52 11.36 11.13 

250 18.06 18.06 18.06 18.05 18.03 17,98 17.92 17.83 17.71 17.57 

350 24.28 24.28 24.28 24.27 24.25 24.21 24.15 24.08 23.98 23.86 

2. = 20 mi/hr 

10 

20 1.92 1.91 1.89 

30 2.23 2.22 2,21 2.14 

50 2.85 2.85 2.83 2.79 2.52 

150 5.96 5.96 5.95 5.93 5.84 5.67 

250 9.07 9.06 9.06 9.05 8.99 8.89 8.73 8.44 

350 12.17 12.17 12.17 12.16 12.12 12.04 11.93 11.76 11.45 

450 15.28 15.28 15.28 15.27 15.23 15.17 15.08 14.94 14.74 

550 18.39 18.39 13.39 18.38 18.35 18.29 18.21 18.09 17.93 17.65 

650 21.50 21.49 21.49 21.19 21.46 21.41 21.33 21.23 21.08 20.85 

750 24.60 24.60 24.60 24.59 24.57 24.52 24.45 24.35 24.21 24.01 

3. v =30 ml/hr 

10 

20 

30 

50 1.92 1.91 1.90 

150 3.99 3.99 3.98 3.95 3.79 

250 6.07 6.06 6.06 6.04 5.95 5.77 

350 8.14 8.14 8.13 8,12 8.06 7.94 7.71 

450 10.21 10.21 10.21 10.19 10.15 10.06 9.91 

550 12.28 12.28 12.28 12.27 12.23 12.16 12.04 11.83 

650 14.35 14.35 14.35 14.34 14.31 14.24 14.14 13.98 

750 16.42 16.42 16.42 16.41 16.38 16.32 16.23 16.09 15.81 

850 18.50 18.50 18.49 18.49 18.46 18.40 18.32 18.19 17.97 

950 20.57 20.57 20.57 20.56 20.53 20.48 20.40 20.28 20.09 

1050 22.64 22.64 22.64 22.63 22.60 22.55 22.48 22.37 22.20 

1150 24.71 24.71 24.71 24.70 24.68 24.63 24.56 24.45 24.28 





URS 702-1 


Table 7A 

SUMMARY OF EXPOSURE DOSE-STANDARD INTENSITY RATIOS 
UP TO H + 36 FOR SELECTED X,Y LOCATIONS 
<w = 1,000 KT, v =20 ml/hr) a 


X (km) D x /I s D(36 )/I s X (km) D l / I s D 2 /I s D ( 36 >/ I s 

Y = 0 km Y = 40 km 

10 0.331 2.54 2.87 150 0.051 1.10 1.15 

20 0.438 2.21 2.65 250 0.033 0.806 0.839 

30 0.446 1.99 2.44 350 0.023 0.620 0.643 

50 0.241 1.70 1.94 450 0.018 0.483 0.501 

150 0.065 1.08 1.15 550 0.014 0.377 0.391 

250 0.038 0.801 0.839 650 0.012 0.289 0.301 

350 0.026 0.617 0.643 750 0.010 0.215 0.225 

450 0.020 0.482 0.502 

550 0.016 0.375 0.391 Y = 60 km 

650 0.013 0.288 0.301 

750 0.011 0.214 0.225 150 0.024 1.12 1.14 

250 0.024 0.813 0.837 

Y = 10 km 350 0.019 0.623 0.642 

450 0.015 0.486 0.501 

10 0.279 2.61 2.89 550 0.012 0.379 0.391 

20 0.410 2.24 2.65 650 0.010 0.291 0.301 

30 0.364 2.01 2.37 750 0.009 0.216 0.225 

50 0.230 1.70 1.93 

150 0.064 1.09 1.15 Y = 80 km 

250 0.038 0.803 0.841 

350 0.026 0.618 0.644 350 0.009 0.632 0.641 

450 0.020 0.482 0.502 450 0.008 0.492 0.500 

550 0.016 0.376 0.392 550 0.007 0.383 0.390 

650 0.013 0.288 0.301 650 0.006 0.294 0.300 

750 0.011 0.214 0.225 750 0.005 0.219 0.224 

Y = 20 km 

30 0.262 2.06 2.32 

50 0.192 1.72 1.91 

150 0.062 1.09 1.15 

250 0.036 0.803 0.839 

350 0.025 0.618 0.643 

450 0.019 0.482 0.501 

550 0.016 0.376 0.392 

650 0.013 0.288 0.301 

750 0.011 0.214 0.225 

a D^ Is the exposure dose to fallout cessation; Y 
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